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Abstract 

In P] John S. Bell proposed how to associate particle trajectories with a lattice 
quantum field theory, yielding what can be regarded as a |^'p-distributed Markov 
process on the appropriate configuration space. A similar process can be defined 
in the continuum, for more or less any regularized quantum field theory; such pro- 
cesses we call Bell-type quantum field theories. We describe methods for explicitly 
constructing these processes. These concern, in addition to the definition of the 
Markov processes, the efficient calculation of jump rates, how to obtain the process 
from the processes corresponding to the free and interaction Hamiltonian alone, 
and how to obtain the free process from the free Hamiltonian or, alternatively, 
from the one-particle process by a construction analogous to "second quantiza- 
tion." As an example, we consider the process for a second quantized Dirac field 
in an external electromagnetic field. 
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1 Introduction 



The aim of this paper is to present methods for constructing Bell-type QFTs. The 
primary variables of Bell-type QFTs are the positions of the particles. Bell suggested a 
dynamical law, governing the motion of the particles, in which the Hamiltonian H and 
the state vector ^ determine certain jump rates [3]. Since these rates are in a sense 
the smallest choice possible, we call them the minimal jump rates. By construction, 
they preserve the distribution. We assume a well-defined Hamiltonian as given; to 
achieve this, it is often necessary to introduce cut-offs. We shall assume this has been 
done where needed. In cases in which one has to choose between several possible position 
observables, for example because of issues related to the Newton- Wigner operator [27| 
1^ . we shall also assume that a choice has been made. 

Bell-type QFTs can also be regarded as extensions of Bohmian mechanics. When 
one tries to incorporate particle creation and annihilation into Bohmian mechanics, one 
is naturally lead to models like the one we presented in ^3]. The quantum equilibrium 
distribution, playing a central role in Bohmian mechanics, then more or less dictates 
that creation of a particle occurs in a stochastic manner — ^just as in Bell's model. 

Bell-type QFTs have in common a good deal of mathematical structure, which we 
will elucidate. The paper is organized as follows. In Section 2 we introduce all the main 
ideas and reasonings; a superficial reading should focus on this section. Some examples 
of Bell-type QFTs are presented in Section 3. (Simple examples of minimal jump rates 
can be found in [IH]-) In Section 4 we describe the construction of a process for the 
free Hamiltonian based on "second quantization." In Section 5 we sketch the concept 
of the "minimal process" associated with a Hamiltonian H . Section 6 concerns some 
properties of Bell-type QFTs that derive from the construction methods developed in 
this paper. In Section 7 we conclude. 



2 Ingredients of Bell- Type Quantum Field Theories 
2.1 Review of Bohmian Mechanics and Equivariance 



Bohmian mechanics [ollTQllzTj is a non-relativistic theory about point particles moving 
in 3-space, according to which the configuration Q = {Qi, . . . ,Qn) evolves according 



toi 



v{Q) , V = nmv 



dt ^ ' ^* ^ ^ ' 

^ = ^t(Q') is the wave function, which evolves according to the Schrodinger equation 

th— = , (2) 



^The masses irik of the particles have been absorbed in the Riemann metric g^i, on configuration 
space R^^, gia,jb = Sij Sab, j — 1 . . . A'^, a, 6 — 1,2,3, and V is the gradient associated with ^aij^i 
i.e., V = {m^^Vq^,...,mj/^VqJ. 
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with 



H 



4a + v 



(3) 



for spinless particles, with A = divV. For particles with spin, ^ takes values in the 
appropriate spin space C'^, V may be matrix valued, and numerator and denominator 
of (P) have to be understood as involving inner products in spin space. The secret 
of the success of Bohmian mechanics in yielding the predictions of standard quantum 
mechanics is the fact that the configuration Qt is | ^-distributed in configuration space 
at all times t, provided that the initial configuration Qq (part of the Cauchy data of the 
theory) is so distributed. This property, called equivariance in [1^], suffices for empirical 
agreement between any quantum theory (such as a QFT) and any version thereof with 
additional (often called "hidden") variables Q, provided the outcomes of all experiments 
are registered or recorded in these variables. That is why equivariance will be our guide 
for obtaining the dynamics of the particles. 

The equivariance of Bohmian mechanics follows immediately from comparing the 
continuity equation for a probability distribution p associated with (^, 



dp 
dt 



- div (/w) , 

with the equation satisfied by |\E'p which follows from (j2I), 



In fact, it follows from (jH]) that 



— Im 

h 



— div 



nim^*(g,t)V^(g,t) 



so, recalling (jT)), one obtains that 



~dt 



-div (1^1 



(4) 



(5) 



(6) 



(7) 



and hence that if pt = l^&tP at some time t then pt = |\l/fp for all times. Equivariance 
is an expression of the compatibility between the Schrodinger evolution for the wave 
function and the law, such as (Q), governing the motion of the actual configuration. In 
[TUj . in which we were concerned only with the Bohmian dynamics (PQ), we spoke of the 
distribution l^'p as being equivariant. Here we wish to find processes for which we have 
equivariance, and we shall therefore speak of equivariant processes and motions. 



2.2 Equivariant Markov Processes 

The study of example QFTs like that of has lead us to the consideration of Markov 
processes as candidates for the equivariant motion of the configuration Q for Hamilto- 
nians H more general than those of the form 
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Consider a Markov process Qt on configuration space. The transition probabilities 
are characterized by the backward generator Lf, a (time-dependent) hnear operator 
acting on functions / on configuration space: 

LJ{q) = -^E{f{Qt^s)\Qt = q) (8) 

where d/ds means the right derivative at s = and E( ■ | ■ ) denotes the conditional 
expectation. Equivalently, the transition probabilities are characterized by the forward 
generator (or, as we shall simply say, generator) , which is also a linear operator but 
acts on (signed) measures on the configuration space. Its defining property is that for 
every process Qt with the given transition probabilities, the distribution pt of Qt evolves 
according to 

is the dual of Lf in the sense that 

J f{q) ^tp{dq) = J Ltfiq) p{dq) . (10) 

We will use both Lt and whichever is more convenient. We will encounter several 
examples of generators in the subsequent sections. 

We can easily extend the notion of equivariance from deterministic to Markov pro- 
cesses. Given the Markov transition probabilities, we say that the |\E'p distribution is 
equivariant if and only if for all times t and t' with t < t', a configuration Qt with dis- 
tribution |\E'tp evolves, according to the transition probabilities, into a configuration Qt' 
with distribution I^E't'p. In this case, we also simply say that the transition probabilities 
are equivariant, without explicitly mentioning Equivariance is equivalent to 

^,i*,r = ^ (11) 

for all t. When (fTT|) holds (for a fixed t) we also say that J/ft is an equivariant generator 
(with respect to \E't and H). Note that this definition of equivariance agrees with the 
previous meaning for deterministic processes. 

We call a Markov process Q equivariant if and only if for every t the distribution pt 
of Qt equals I^E'tp. For this to be the case, equivariant transition probabilities are nec- 
essary but not sufficient. (While for a Markov process Q to have equivariant transition 
probabilities amounts to the property that if pj = |\E'tp for one time t, where pt denotes 
the distribution of Qt, then p^/ = l^^t'l"^ for every t' > t, according to our definition of 
an equivariant Markov process, in fact pt = l^&tP for all t.) However, for equivariant 
transition probabilities there exists a unique equivariant Markov process. 

The crucial idea for our construction of an equivariant Markov process is to note 
that is completely general, and to find a generator J^t such that the right hand side 
of (0) can be read as the action of ^ on p = |\E'p, 

|lm**/7^ = ^l^p. (12) 
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We shall implement this idea beginning in Section after a review of jump processes 
and some general considerations. But first we shall illustrate the idea with the familiar 
case of Bohmian mechanics. 

For H of the form Q, we have ^ and hence that 

2 / n \I^*V^\ 

-lmm*Hm = -div (;iIm^*V^) = - div M^p^iIm^^J . (13) 

Since the generator of the (deterministic) Markov process corresponding to the dynam- 
ical system dQ/dt = v{Q) given by a velocity vector field v is 

^p = -divipv), (14) 

we may recognize the last term of (jl3|) as =Sf |\E'p with =Sf the generator of the determin- 
istic process defined by (Q). Thus, as is well known, Bohmian mechanics arises as the 
natural equivariant process on configuration space associated with H and 

To be sure, Bohmian mechanics is not the only solution of (fT^ for H given by 
Q. Among the alternatives are Nelson's stochastic mechanics and other velocity 
formulas [12j. However, Bohmian mechanics is the most natural choice, the one most 
likely to be relevant to physics. It is, in fact, the canonical choice, in the sense of minimal 
process which we shall explain in Section 15.31 

2.3 Equivariant Jump Processes 

Let Q denote the configuration space of the process, whatever sort of space that may 
be (vector space, lattice, manifold, etc.); mathematically speaking, we need that Q be a 
measurable space. A (pure) jump process is a Markov process on Q for which the only 
motion that occurs is via jumps. Given that Qt = q, the probability for a jump to q', 
i.e., into the infinitesimal volume dq' about q', by time t + dt is (Tt{dq'\q) dt, where a is 
called the jump rate. In this notation, cr is a finite measure in the first variable; a{B\q) 
is the rate (the probability per unit time) of jumping to somewhere in the set B C 
given that the present location is q. The overall jump rate is (j{Q\q). 

It is often the case that Q is equipped with a distinguished measure, which we shall 
denote by dq or dq', slightly abusing notation. For example, if Q = M.'^, dq may be the 
Lebesgue measure, or if Q is a Riemannian manifold, dq may be the Riemannian volume 
element. When cr( ■ |g) is absolutely continuous relative to the distinguished measure, 
we also write o"(g'|g) dq' instead of a{dq'\q). Similarly, we sometimes use the letter p for 
denoting a measure and sometimes the density of a measure, p{dq) = p{q) dq. 

A jump first occurs when a random waiting time T has elapsed, after the time to at 
which the process was started or at which the most recent previous jump has occurred. 
For purposes of simulating or constructing the process, the destination q' can be chosen 
at the time of jumping, to + T, with probability distribution (Ttg+T{Q\q)~^ <^to+T{ ■ \q)- 
In case the overall jump rate is time-independent, T is exponentially distributed with 
mean cr(Q|g)~^. When the rates are time- dependent — as they will typically be in what 
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follows — the waiting time remains such that 

rto+T 

/ cTtiQ\q)dt 
J to 

is exponentially distributed with mean 1, i.e., T becomes exponential after a suitable 
(time-dependent) rescaling of time. For more details about jump processes, see jH]. 
The generator of a pure jump process can be expressed in terms of the rates: 

^ap{dq) = J ((y{dq\q')p{dq') - a{dq'\q)p{dq)^ , (15) 

q'eQ 

a "balance" or "master" equation expressing dp/dt as the gain due to jumps to dq minus 
the loss due to jumps away from q. 

We shall say that jump rates a are equivariant if ^„ is an equivariant generator. It 
is one of our goals in this paper to describe a general scheme for obtaining equivariant 
jump rates. In Sections 12.61 and ITTI we will explain how this leads us to formula (1^. 



2.4 Process Additivity 

The Hamiltonian of a QFT usually comes as a sum, such as 

H = Ho + Hi (16) 

with Ho the free Hamiltonian and Hj the interaction Hamiltonian. If several particle 
species are involved, Hq is itself a sum containing one free Hamiltonian for each species. 
The left hand side of (jl2p . which should govern our choice of the generator, is then also 
a sum, 

-Im^*iyo^ + -Im$*if/^ = ^|^P. (17) 
n h 

This opens the possibility of finding a generator =Sf by setting =Sf = + =Sf/, provided 
we have generators and corresponding to Hq and Hj in the sense that 

^ol^P (18a) 
^/l^l^ (18b) 

This feature of ()12|1 we call process additivity; it is based on the fact that the left hand 
side of fll2|l is linear in H. Note that the backward generator of the process with forward 
generator =Sfo + =Sf/ is Lq + Lj; thus forward and backward generators lead to the same 
notion of process additivity, and to the same process corresponding to Hq + Hj. In 
many cases, as will be elaborated in Section Hq is based on an operator known from 
quantum mechanics (e.g., the Dirac operator), in such a way that =Sfo can be obtained 
from the appropriate Bohmian law of motion. In Section 12.61 we will explain how =2/ 
can usually be taken as the generator of a jump process. 



n 

- Im ^!*Hi^! = 

h 
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Our proposal is to take seriously the process generated by =Sf = ££q + and regard 
it as the process naturally associated with H. The bottom line is that process additivity 
provides a method of constructing a Bell-type theory. 

Obviously, the mathematical observation of process additivity (that sums of gener- 
ators define an equivariant process associated with sums of Hamiltonians) applies not 
only to the splitting of H into a free and an interaction contribution, but to every case 
where if is a sum. And it seems that process additivity provides a physically very rea- 
sonable process in every case where H is naturally a sum, in fact the most reasonable 
process: the one that should be considered the Bell-type process, defining the Bell-type 
theory. 

2.5 What Added Processes May Look Like 

To get some feeling for what addition of generators, ^ = =Sfi + J^2, means for the cor- 
responding processes, we consider some examples. First consider two deterministic pro- 
cesses (on the same configuration space), having generators of the form =Sfp = — div {pv). 
To add the generators obviously means to add the velocity vector fields, f = f i + f2, so 
the resulting velocity is a superposition of two contributions. 

Next consider a pure jump process. Since, according to (fT3j) . the generator ^ is 
linear in a, adding generators means adding rates, a = (Xi + cr2- This is equivalent to 
saying there are two kinds of jumps: if the present location is g G Q, with probability 
CiiQlq) dt the process performs a jump of the first type within the next dt time units, and 
with probability cr2{Q\q) dt a jump of the second type. That does not mean, however, 
that one can decide from a given realization of the process which jump was of which 
type. 

Next suppose we add the generators of a deterministic and a jump process, 

ifp(g) = - div (pf )(g) + J (cT(g|g') p{q') - a{q'\q) p{q)^dq . (19) 

g'eS 

This process moves with velocity v{q) until it jumps to g', where it continues moving, 
with velocity v{q'). The jump rate may vary with time in two ways: first because a 
may be time-dependent, second because a may be position-dependent and Qt moves 
with velocity v. One can easily understand ()19p in terms of gain or loss of probability 
density due to motion and jumps. So this process is piecewise deterministic: although 
the temporal length of the pieces (the intervals between two subsequent jumps) and the 
starting points (the jump destinations) are random, given this data the trajectory is 
determined. 

The generator of the Wiener process in is the Laplacian, and to add to it the 
generator of a deterministic process means to introduce a drift. Note that this is different 
from adding, in W^, a Wiener process to a solution of the deterministic process. In spaces 
like W^, where it so happens that one is allowed to add locations, there is a danger of 
confusing addition of generators with addition of realizations. Whenever we speak of 
adding processes, it means we add generators. 
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To add generators of a diffusion and a pure jump process yields what is often called a 
jump diffusion process, one making jumps with time- and position-dependent rates and 
following a diffusion path in between. Diffusion processes, however, will play almost no 
role in this paper. 



2.6 Integral Operators Correspond to Jump Processes 

We now address the interaction part Hj of the Hamiltonian (fTHI) . In QFTs with cutoffs 
it is usually the case that Hi is an integral operator. For that reason, we shall in this 
work focus on integral operators for Hi. We now point out why the naturally associated 
process is a pure jump process. For short, we will write H rather than Hi in this and 
the subsequent section. For the time being, think of Q as and of wave functions as 
complex valued. 

What characterizes jump processes versus continuous processes is that some amount 
of probability that vanishes at g G Q can reappear in an entirely different region of 
configuration space, say at q' G Q. This is manifest in the equation for dp/dt, (fT3j): 
the first term in the integrand is the probability increase due to arriving jumps, the 
second the decrease due to departing jumps, and the integration over q' reflects that q' 
can be anywhere in Q. This suggests that Hamiltonians for which the expression (0) 
for d\^\'^ /dt is naturally an integral over dq' correspond to pure jump processes. So 
when is the left hand side of ()12|) an integral over dq'l When H is an integral operator, 
i.e., when {q\H\q') is not merely a formal symbol, but represents an integral kernel that 
exists as a function or a measure and satisfies 

{Hm){q) = j dq'{q\H\q')<il{q'). (20) 



In this case, we should choose the jump rates in such a way that, when p = |\E' 



2 



aiq\q') p(g') - a(g'|g) p{q) = ^ Im ^*iq) {q\H\q') ^{q') , (21) 



and this suggests, since jump rates must be nonnegative (and the right hand side of (j2T 
is anti-symmetric), that 

criqW) Piq') = [^Imvl/*(g) {q\H\q') ^{q') 

(where x"*" denotes the positive part of x G M, that is, x"*" is equal to x for x > and is 
zero otherwise), or 



, I [{2/h)lm^*{q) {q\H\q')^{q'), 

"(^'^ ) = ^^mw) • ^''^ 

These rates are an instance of what we call the minimal jump rates associated with H 
(and \E'). The name comes from the fact that they are actually the minimal possible 
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values given (HH), as is expressed by the inequality ()114|) and will be explained in detail 
in Section 15.21 Minimality entails that at any time t, one of the transitions qi —>■ q2 or 
?2 Qi is forbidden. We will call the process defined by the minimal jump rates the 
minimal jump process (associated with H). 

In contrast to jump processes, continuous motion, as in Bohmian mechanics, cor- 
responds to such Hamiltonians that the formal matrix elements {q\H\q') are nonzero 
only infinitesimally close to the diagonal, and in particular to differential opera- 
tors like the Schrodinger Hamiltonian which has matrix elements of the type 
6"{q — q') + V{q)5{q — q'). We can summarize the situation, as a rule of thumb, by 
the following table: 



A contribution to H that is a . . . 


corresponds to . . . 


integral operator 
differential operator 
multiplication operator 


jumps 

deterministic continuous motion 
no motion = 0) 



The minimal jump rates as given by ()22|) have some nice features. The possible 
jumps for this process correspond to the nonvanishing matrix elements of H (though, 
depending on the state ^f, even some of the jump rates corresponding to nonvanishing 
matrix elements of H might happen to vanish). Moreover, in their dependence on the 
state the jump rates a depend only "locally" upon the jump rate for a given jump 
q' q depends only on the values '^{q') and \E'(g) corresponding to the configurations 
linked by that jump. Discretizing M? to a lattice el?, one can obtain Bohmian mechanics 
as a limit e — > of minimal jump processes [331 whereas greater-than-minimal jump 
rates lead to Nelson's stochastic mechanics [2E! and similar diffusions, such as ()117|1 : 
see I311I22!- If the Schrodinger operator (jH} is approximated in other ways by operators 
corresponding to jump processes, e.g., by H;, = e~^^ He~^^ , the minimal jump processes 
presumably also converge to Bohmian mechanics. 

We have reason to believe that there are lots of self-adjoint operators which do not 
correspond to any stochastic process that can be regarded as defined, in any reasonable 
sense, by fE^.^ But such operators seem never to occur in QFT. (The Klein-Gordon 
operator ^/m^c^ — K^c^A does seem to have a process, but it requires a more detailed 
discussion which will be provided in a forthcoming work [TH].) 

2.7 Minimal Jump Rates 

The reasoning of the previous section applies to a far more general setting than just con- 
sidered: to arbitrary configuration spaces Q and "generalized observables" — POVMs — 
defining, for our purposes, what the "position representation" is. We now present this 
more general reasoning, which leads to one of the main formulas of this paper, ()29|) . 

^Consider, for example, H = pcosp where p is the one-dimensional momentum operator —ihd/dq. 
Its formal kernel {q\H\q') is the distribution —^S'{q — g' — 1) — ^S'{q — q' + 1), for which would not 
have a meaning. From a sequence of smooth functions converging to this distribution, one can obtain 
a sequence of jump processes with rates (I22|) : the jumps occur very frequently, and are by amounts of 
approximately ±1. A limiting process, however, does not exist. 
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The process we construct relies on the following ingredients from QFT: 

1. A Hilbert space with scalar product (^E'l^). 

2. A unitary one-parameter group \Jt in ^ with Hamiltonian if, 

so that in the Schrodinger picture the state evolves according to 

Ut could be part of a representation of the Poincare group. 

3. A positive-operator-valued measure (POVM) P{dq) on Q acting on J^, so that 
the probability that the system in the state \& is localized in dq at time t is 

Ftidq) = {^t\P{dq)\^t) . (24) 

Mathematically, a POVM P on Q is a countably additive set function ("measure"), 
defined on measurable subsets of Q, with values in the positive (bounded self-adjoint) 
operators on (a Hilbert space) Jif, such that P{Q) is the identity operator.^ Physically, 
for our purposes, P{-) represents the (generalized) position observable, with values in 
Q. The notion of POVM generalizes the more familiar situation of observables given 
by a set of commuting self-adjoint operators, corresponding, by means of the spectral 
theorem, to a projection- valued measure (PVM): the case where the positive operators 
are projection operators. A typical example is the single Dirac particle: the position 
operators on L'^{M.^,C'^) induce there a natural PVM Po{-)'- for any Borel set B C 
M^, Po{B) is the projection to the subspace of functions that vanish outside B, or, 
equivalently, Po{B)'^{q) = ls(?) ^(q') with 1b the indicator function of the set B. 
Thus, (\l/|Po('^Q')|^) = \'^{(l)\'^dq. When one considers as Hilbert space Jif only the 
subspace of positive energy states, however, the localization probability is given by 
P( . ) = p+Po( . )/ with P+ : L2(M^ C^) ^ the projection and / : ^ ^ L\R^, C^) 
the inclusion mapping. Since P+ does not commute with most of the operators Pq{B), 
P( ■ ) is no longer a PVM but a genuine POVM^ and consequently does not correspond to 
any position operator — although it remains true (for \E' in the positive energy subspace) 
that (\E'|P((ig)|^) = |\E'(g)p(ig. That is why in QFT, the position observable is indeed 
more often a POVM than a PVM. POVMs are also relevant to photons fP, In one 
approach, the photon wave function \E' : — is subject to the constraint condition 
V ■ \E' = di'^i + d2'^2 + ds^s = 0. Thus, the physical Hilbert space is the (closure 

■^The countable additivity is to be understood as in the sense of the weak operator topology. This 
in fact implies that countable additivity also holds in the strong topology. 

''This situation is indeed more general than it may seem. By a theorem of Naimark ^2 P- 142], 
every POVM P{ ■ ) acting on is of the form P( • ) = P+Po( • )P+ where Pq is a PVM on a larger 
Hilbert space, and P+ the projection to J^. 
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of the) subspace of L^(M^,C^) defined by this constraint, and the natural PVM on 
L^(M^,C^) gives rise, by projection, to a POVM on . So much for POVMs. Let us 
get back to the construction of a jump process. 

The goal is to specify equivariant jump rates o = cr*'^''^, i.e., such rates that 

r/P 

XP=^. (25) 

To this end, one may take the following steps: 

1. Note that 

^^p^ = llm{%\P{dq)H\%). (26) 

2. Insert the resolution of the identity / = J P{dq') and obtain 

g'eQ 



dt 

g'eQ 



{dq^dq'), (27) 



where 

Udq,dq') = "^Im {^t\Pidq)HP{dq')\^t) • (28) 

3. Observe that J is anti-symmetric, S{dq',dq) = —S{dq,dq'). Thus, since x = — 

(-^)^ 

I{dq,dq') = [{2/h) lm{^\P{dq)HP{dq')\^)f - [{2/h) lm{^\P{dq')HP{dq)\^)f . 

4. Multiply and divide both terms by P( • ), obtaining that 

g'eQ g'eS 

[{2/h)lm{^\P{dq')HP{dqm] + ^^^^^\ 
{^\P{dq)\^) "^^'^V ■ 

5. By comparison with (fT3j) . recognize the right hand side of the above equation as 
=Sfo-P, with =Sfo- the generator of a Markov jump process with jump rates 

^(rJn\n'\ - iiV ^) ^^^^ {^\P idq)H P idq')\^)]^ 

{^\P{dq')\^) ' ^ 

which we call the minimal jump rates. 
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Mathematically, the right hand side of this formula as a function of q' must be understood 
as a density (Radon-Nikodym derivative) of one measure relative to another.^ The plus 
symbol denotes the positive part of a signed measure; it can also be understood as 
applying the plus function, x"*" = max(x, 0), to the density, if it exists, of the numerator. 

To sum up, we have argued that with H and ^ is naturally associated a Markov 
jump process Qt whose marginal distributions coincide at all times by construction with 
the quantum probability measure, Pt{-) = '^ti')-, so that Qt is an equivariant Markov 
process. 

In Section 4 of we establish precise conditions on iJ, P, and \E' under which the 
jump rates are well-defined and finite P-almost everywhere, and prove that in this 
case the rates are equivariant, as suggested by the steps 1-5 above. It is perhaps worth 
remarking at this point that any H can be approximated by Hamiltonians (namely 
Hilbert-Schmidt operators) for which the rates (j29|) are always (for all ^) well-defined 
and equivariant [T3] . Concerning this, see also the end of Section 15.31 



2.8 Process Associated with the Free Hamiltonian 

We now address the free Hamiltonian Hq of a QFT. We describe the process naturally 
associated with Hq, when this is the second quantized Schrodinger or Dirac operator. 
We will treat more general free Hamiltonians in the next section. We shall consider here 



^Quite aside from the previous discussion, it is perhaps worth noting that there are not so many 
expressions in H, P, and \1/ that would meet the formal criteria for being a candidate for the jump 
rate. Since the only connection between abstract Hilbert space and configuration space is by P, which 
leads to measures on Q, the only way to obtain a function on Q is to form a Radon-Nikodym quotient 
of two measures, (j{q') = A{dq') / B{dq'). Since a must be a measure- valued function, the numerator 
should be a bi-measure (a measure in each of two variables). The simplest measure one can form from 
i/,P, and is (*|P(dg)|^'); the simplest bi-measures are (*|iJ"iP((ig)7?"2P(dg')-ff"^ I*)- Jump rates 
must have dimension 1/time, and the only object at hand having this dimension is H/h. Thus, H can 
appear only once in the numerator. The expressions {'^\HP{dq)P{dq')Y^) and {^\P{dq)P{dq')H\^) 
are no good because for PVMs P they are concentrated on the diagonal of Q x Q and hence do not 
lead to nontrivial jumps. Let us write /i for the measure-valued function we have arrived at: 

1 {^P{dq)HP{dq'm 

^('^'^''^^ = fi mpidq'm ■ 

This provides complex measures, whereas cr( • \q') must be a positive real measure. There are not many 
ways of forming a positive real measure from a complex one, the essential ones being 

1^1, |Re^|, |Im^|, (Re^)+, (Re^)", (Im^)+, (Im/^)" 

times a numerical constant A > 0. One could of course form additional expressions at the price of 
higher complexity. 

This has gotten us already pretty close to the minimal rates (|29ll . which correspond to cr = 2(Im/i)+. 
To proceed further, we might demand the absence of unnecessary jumps; that means that at any time, 
either the jump qi — > q2 or q2 — > q\ is forbidden; this leaves only A(Im/x)^. Moreover, 2(Im/i)+ is 
the only expression in the list that has Bohmian mechanics as a limiting case or implies equivariance. 
Furthermore it corresponds to the natural guess (|118|l for a backward generator, discussed in Section 
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only Hamiltonians for one type of particle. 

We first define the configuration space Q. Let us write Q^^^ ("one-particle configu- 
ration space") for physical space; this is typically, but not necessarily, R^. The space 
Q in which the "free process" takes place is the configuration space for a variable num- 
ber of identical particles; we call it TQ.^^\ It can be defined as the space of all finite 
subsets- with-multiplicities of QP''^. A set- with- multiplicities consists of a set and, for 
each element x of the set, a positive integer, called the multiplicity of x. The number 
of particles in a configuration q is the sum of its multiplicities, Such configura- 
tions describe several identical particles, some of which may be located at the same 
position in space. Equivalently, one could say that FQ*^^) is the set of all mappings 
n : 

^ Nu {0} (meaning the number of particles at a given location) such that 

^ n{q) < oo. 

qeQW 

Another equivalent definition is the set of all finite nonnegative measures n( ■ ) on Q^^^ 
that assume only integer values; the meaning of n{R) is the number of particles in the 
region R of physical space. Finally, one can define 

oo 

FQ(^) = IJ where Q^") = (Q(^))"/permutations. 

n=0 

A related space, for which we write F^Q^^^ is the space of all finite subsets of Q^^^; 
it is contained in TQ^^\ after obvious identifications. In fact, T^Q^^^ = TQ'^^^ \ A, where 
A is the set of coincidence configurations, i.e., those having two or more particles at the 
same position. T^Q^^^ is the union of the spaces Q^"* for n = 0, 1, 2, . . ., where Q^T'' is 
the space of subsets of Q^^^ with n elements. 

For Q^^^ = W^, the n-particle sector Q^"* is a manifold of dimension nd (see jTB] for 
a discussion of Bohmian mechanics on this manifold). If d > 2, the set A of coincidence 
configurations has codimension > 2 and thus can usually be ignored. We can then 
replace FM'^ by the somewhat simpler space F^M''. 

The position POVM P^^^ on Q^^^ (acting on the one-particle Hilbert space) naturally 
leads to a POVM we call TP^^^ on Q = TQ^^\ acting on Fock space (see Section ll.2.3l for 
the definition).^ Since a configuration from F(M^) defines the number of particles and 
their positions, the name "position observable" for P = TP'^^^ stretches the meaning of 
"position" somewhat: it now also encompasses the number of particles. 

We now give a description of the free process associated with the second- quantized 
Schrodinger operator; it arises from Bohmian mechanics. Fock space Jt^ = ^ is a direct 
sum 

oo 

^ = 0^("), (30) 

n=0 

^The coincidence configurations form a null set, rp(^)(A) = 0, when Q'^^^ is a continumir, or, more 
precisely, when P^'^-' is nonatomic as a measure. 
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where ^("^ is the ?7,-particle Hilbert space. ^("^ is the subspace of symmetric (for 
bosons) or anti-symmetric (for fermions) functions in L^(M^", (^C^'*"''^)'^") for spin-s par- 
ticles. Thus, \I/ G ^ can be decomposed into a sequence \1/ = . . . , . . 
the n-th member being an ri-particle wave function, the wave function represent- 
ing the n-particle sector of the quantum state vector. The obvious way to obtain a 
process on Q = FM^ is to let the configuration Q{t), containing = 4i^Q{t) particles, 
move according to the A^-particle version of Bohm's law (0), guided by \['(^).'^ This is 
indeed an equivariant process since Hq has a block diagonal form with respect to the 
decomposition (jHOJ, 

oo 
n=0 

and -f^o"'* is just a Schrodinger operator for n noninteracting particles, for which, as we 
already know, Bohmian mechanics is equivariant. We used a very similar process in J3] 
(the only difference being that particles were numbered in [T^). 

Similarly, if Hq is the second quantized Dirac operator, we let a configuration Q with 
particles move according to the usual A^-particle Bohm-Dirac law [7, p. 274] 

dQ^ vl/*(g) aN^{Q) . . 

dt ^*{Q)^{Q) ^ ^ 

where c denotes the speed of light and = {ol'^^\ • • • , ol^^^) with a'^'^^ acting on the 
spin index of the /c-th particle. 



2.9 Other Approaches to the Free Process 

We will give below a general velocity formula, applicable to a wider class of free Hamil- 
tonians. Alternatively, we can provide a free process for any Hq if we are given an 
equivariant process for the one-particle Hamiltonian H'^^\ This is based on the par- 
ticular mathematical structure of Hq, which can be expressed by saying it arises from 
a one-particle Hamiltonian H^^^ by applying a "second quantization functor F" |29j . 
That is, there is an algorithm (in a bosonic or fermionic version) for forming, from 
a one-particle Hilbert space J^^^^ and a one-particle Hamiltonian H^^\ a Fock space 
^ = TJif^^^ and free Hamiltonian Hq = TH^^\ And parallel to this "second quanti- 
zation" algorithm, there is an algorithm for the canonical construction, from a given 
equivariant one-particle Markov process Q[^\ of a process we call FQI^^ that takes place 
in Q = TQ^^^ and is equivariant with respect to Hq. This algorithm may be called the 
"second quantization" of a Markov process. 

The algorithm is described in Section 14.21 What the algorithm does is essentially to 
construct an n-particle version of gl^^ for every n, and finally combine these by means 
of a random particle number A^ = N{t) = i^Q{t) which is constant under the free 

^As defined, configurations are unordered, whereas we fiave written Bolim's law ^ for ordered 
configurations. Tlianks to tlie (anti-)symmetry of the wave function, however, all orderings will lead to 
the same particle motion. For more about such considerations, see our forthcoming work jl3|. 



15 



process, parallel to the fact that the particle number operator is conserved by Hq. We 
note further that the process TQ[^^ is deterministic if Qf^ is. If we take the one-particle 
process to be Bohmian mechanics or the Bohm-Dirac motion, the algorithm reproduces 
the processes described in the previous section. 

The algorithm leaves us with the task of finding a suitable one-particle law, which 
we do not address in this paper. For some Hamiltonians, such as the Dirac operator, 
this is immediate, for others it is rather nontrivial, or even unsolved. The Klein-Gordon 
operator y/m'^c'^ — K^c^A will be discussed in forthcoming work [TH] , and for a study of 
photons see [2H]- 

When Ho is made of differential operators of up to second order (which includes of 
course the Schrodinger and Dirac operators), there is another way to characterize the 
process associated with Hq, a way which allows a particularly succinct description of the 
process and a particularly direct derivation and construction. In fact, we give a formula 
for its backward generator Lq, or alternatively the velocity (or the forward generator 
=Sfo), in terms of Hq, P, and \E'. 

We begin by defining, for any H, P, and an operator L acting on functions 
/ : Q — i> R, which may or may not be the backward generator of a process, by 

where [ , 1 means the commutator, 



f= J f{q)P{dq), (33) 

q&Q 

and L is the "generator" of the (Heisenberg) time evolution of the operator /, 



dr 



T=0 



mj]. (34) 



(If P is a PVM, then / = /(g), where q is the configuration operator.) (jH^ could be 
guessed in the following way: since Lf is in a certain sense, see (jH)), the time derivative 
of /, it might be expected to be related to Lf, which is in a certain sense, see dSH), the 
time derivative of /. As a way of turning the operator Lf into a function Lf{q), the 
middle term in (jH^ is an obvious possibility. Note that this way of arriving at (jH^ does 
not make use of equivariance; for another way that does, see Section 15.11 
The formula for the forward generator equivalent to (j32j) reads 

^p{dq) = Re{^\%j:[H,P{dq)m, (35) 

as follows from (fTUj) . 

Whenever L is indeed a backward generator, we call it the minimal free (backward) 
generator associated with '^,H, and P. (The name is based on the concept of minimal 
process as explained in Section 1^31 ) Then the corresponding process is equivariant (see 
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Section EH}. This is the case if (and, there is reason to expect, only if) P is a PVM 
and if is a differential operator of up to second order in the position representation, 
in which P is diagonal In that case, the process is deterministic, and the backward 
generator has the form L = v -V where v is the velocity vector field; thus, ()32j) directly 
specifies the velocity, in the form of a first-order differential operator v ■ V. In case H 
is the A^-particle Schrodinger operator with or without spin, ()32j) yields the Bohmian 
velocity (0), and if H is the Dirac operator, the Bohm-Dirac velocity (j31|) . To sum up, 
in some cases definition leads to just the right backward generator. 

To return to our starting point: if the one-particle generator Sf^^^ arises from the 
one-particle Hamiltonian H^^^ by (jH^j) . then (jH^j) also holds between the free generator 
=5fo = T^^^^ and the free Hamiltonian Hq = TH^^\ (See Section f5.ll for details.) In 
other words, (j32j) is compatible with the "second quantization" algorithm. Thus, in 
relevant cases (jH^ allows a direct definition of the free process in terms of Hq, just as 
(Uni) directly defines, in terms of Hj, the jump rates. 

A relevant point is that the "second quantization" of a differential operator is again 
a differential operator, in a suitable sense, and has the same order. Note also that ()32|). 
when applied to the second quantized Schrodinger or Dirac Hamiltonian, defines the 
same vector field on r(M^) as described in the previous section. 

2.10 Bell-Type QFT 

We briefly summarize what we have obtained. A Bell-type QFT is about particles 
moving in physical 3-space; their number and positions are represented by a point Qt in 
configuration space Q. Provided physical space is M^, Q is usually FM^ or a Cartesian 
product of several such spaces, each factor representing a different particle species. Qt 
follows a Markov process in Q, which is governed by a state vector \E' in a suitable 
Hilbert space J^. M' is related to Q, by means of a PVM or POVM P. \E' undergoes a 
unitary evolution with Hamiltonian H . The process Qt usually consists of deterministic 
continuous trajectories interrupted by stochastic jumps; more generally, it arises by 
process additivity (i.e., by adding generators) from a free process associated with and 
a jump process associated with Hj. The jump rates are given by ()29|) for H = Hj. The 
free process arises from Bohmian mechanics, or a suitable analogue, by a construction 
that can be formalized as the "second quantization" of a one-particle Markov process; 
when appropriate, it is defined directly by (j^ . The process Qt is equivariant, i.e., 
(^t|P(dg)|^t) distributed. 

Examples of Bell-type QFTs can be found in [3 [Ej and in Section El It is our 
contention that, essentially, there is a unique Bell- type version of every regularized 
QFT. We have to postpone, however, the discussion of operators of the Klein-Gordon 
type. We also have to assume that the QFT provides us with the POVM P{-)', this 
is related to an ongoing discussion in the literature j2Zl 1211 12^1 concerning the right 
position operator. 
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2.11 More on Identical Particles 



The n-particlc sector of the configuration space (without coincidence configurations) 
of identical particles r^(R"^) is the manifold of n-point subsets of M^; let Q be this 
manifold. The most common way of describing the quantum state of n fermions is by 
an anti-symmetric (square- integrable) wave function on Q := R^"; let be the 
space of such functions. Whereas for bosons ^ could be viewed as a function on Q, for 
fermions is not a function on Q. 

Nonetheless, the configuration observable still corresponds to a PVM P on Q: for 
B C Q, we set P{B)^{q^, ■■■,qn) = ■ ■ • , 9n) if {^i, • • ■ , q„} G 5 and zero oth- 

erwise. In other words, P{B) is multiplication by the indicator function of t:~^{B) 
where tt is the obvious projection mapping Q \ A ^ Q, with A the set of coincidence 
configurations. 

To obtain other useful expressions for this PVM, we introduce the formal kcts \q) 
for g G Q (to be treated like elements of L^(Q)), the anti-symmetrization operator S 
(i.e., the projection LF'{Q) ^), the normalized anti-symmetrizer^ s — Vn^S, and 
the formal kets \sq) :— s\q) (to be treated like elements of J^f). The \q) and are 
normalized in the sense that 

{q\q') = 5{q - q') and {sq\sq') = (-1)^?«,5') s{q - g'), 

where q = 7i{q), q' = 7i{q'), g{q,q') is the permutation that carries q into q' given that 
q — q'i and (—1)^ is the sign of the permutation q. Now we can write 

P{dq)= \ot){q\dq = n\S\q){q\dq=\sq){sq\dq, (36) 

?e7r-i(q) 

where the sum is over the n! ways of numbering the n points in g; the last two terms 
actually do not depend on the choice of ^ e 7r~^(gf), the numbering of q. 
The probability distribution arising from this PVM is 

^{dq)= \-<^{q)\'^dq = n\\-<^{q)\'^dq=\{sq\-<^)\^dq (37) 

with arbitrary q G 7r^^(g). 

There is a way of viewing fermion wave functions as being defined on Q, rather than 
M^", by regarding them as cross- sections of a particular 1-dimensional vector bundle 
over Q. To this end, define an n!-dimensional vector bundle E by 

E,:= C. (38) 

qe-K-^(q) 

^The name means this: since 5 is a projection, S'^ is usually not a unit vector when is. Whenever 
^ G L^{Q) is supported by a fundamental domain of the permutation group, i.e., by a set O C Q on 
which (the restriction of) tt is a bijection to Q, the norm of is l/\/n\, so that s"^ is again a unit 
vector. 
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Every function \1/ : M.^^ C naturally gives rise to a cross-section $ of E, defined by 

$(g):= ^(q). (39) 

ge7r-l{g) 

The anti-symmetric functions form a 1-dimensional subbundle of E (see also 13. for a 
discussion of this bundle) . 

3 Application to Simple Models 

In this section, we point out how the jump rates of the model in are contained in 
()29|) and present a full-fledged Bell-type QFT for the second-quantized Dirac equation 
in an external electromagnetic field. 

Further cut-off QFTs that may provide interesting examples of Bell-type QFTs, 
worth a detailed discussion in a future work ^7j, are the scalar self-interacting field 
(e.g., $^), QED, and other gauge field theories. We have to postpone the treatment 
of these theories because they require discussions lying outside the scope of this paper, 
in particular a discussion of the position representation of photon wave functions in 
QED, and, concerning of the appropriate probability current for the Klein-Gordon 
equation. 

3.1 A Simple QFT 

We presented a simple example of a Bell-type QFT in and we will now briefly point 
to the aspects of this model that are relevant here. The model is based on one of the 
simplest possible QFTs [221 P- 339]. 

The relevant configuration space Q for a QFT (with a single particle species) is 
the configuration space of a variable number of identical particles in M^, which is the 
set r(M^), or, ignoring the coincidence configurations (as they are exceptions), the set 
r^(M^) of all finite subsets of M^. The n-particle sector of this is a manifold of dimension 
3n; this configuration space is thus a union of (disjoint) manifolds of different dimensions. 
The relevant configuration space for a theory with several particle species is the Cartesian 
product of several copies of r^(M^). In the model of ^14j, there are two particle species, 
a fermion and a boson, and thus the configuration space is 

Q = T^{R^) X T^{M.^). (40) 

We will denote configurations by g = (x, y) with x the configuration of the fermions and 
y the configuration of the bosons. 

For simplicity, we replaced in J3] the sectors of (M^) x F^ (M^) , which are manifolds, 
by vector spaces of the same dimension (by artificially numbering the particles), and 
obtained the union 

oo oo 

Q = |J(M3)n X [j (M^)™, (41) 
n=0 m=0 
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with n the number of fermions and m the number of bosons. Here, however, we will use 
(jiUj) as the configuration space, since we have already discussed the space r^(R'^). In 
comparison with (j4H) . this amounts to (merely) ignoring the numbering of the particles. 

Jif is the tensor product of a fermion Fock space and a boson Fock space, and thus the 
subspace of wave functions in ( Q) that are anti-symmetric in the fermion coordinates 
and symmetric in the boson coordinates. Let S denote the appropriate symmetrization 
operator, i.e., the projection operator L'^{Q) J^, and s the normalized symmetrizer 

s'${xi,...,Xn,yi,...,y^) = Vn\ m\ S^{xi, . . . , a;„, y^, . . . , t/^), (42) 

i.e., s = a/ A^! M! 5* with and M the fermion and boson number operators, which 
commute with S and with each other. As in Section r2.1H we denote by vr the projection 
mapping Q \ A ^ Q, ti{xi, . . . , a;„, y^, . . . , y„) = {{xi, . . . , x^}, {2/1, . . . , ?/„}). The 
configuration PVM P{B) on Q is multiplication by 17^-1(5), which can be understood as 
acting on J^, though it is defined on L^(Q), since it is permutation invariant and thus 
maps to itself. We utilize again the formal kets \q) where g G Q \ A is a numbered 
configuration, for which we also write q = {x, y) = {xi, . . . , Xn, Hi, . . . , l/m)- We also use 
the symmetrized and normalized kets \sq) = s\q). As in ()36|) . we can write 

P{dq) = \q){q\dq = n\m\ S\q){q\dq = \sq){sq\dq (43) 

ge7r-i((j) 

with arbitrary q G 7r~^(g). For the probability distribution, we thus have, as in ()37|1 . 

F{dq)= |*(g)|2c/g = n!m! |*(g)|^rfg = |(sg|^)|2c/g (44) 

(?e7r-i{g) 

with arbitrary q G 7r~^(g). 

The free Hamiltonian is the second quantized Schrodinger operator (with zero poten- 
tial), associated with the free process described in Section 17!^ The interaction Hamil- 
tonian is defined by 

Hj = j d^xi,\x){al{x) + a^{x))i,{x) (45) 

with il)'^{x) the creation operators (in position representation), acting on the fermion 
Fock space, and a^(a;) the creation operators (in position representation), acting on the 
boson Fock space, regularized through convolution with an function : — M. Hj 
has a kernel; we will now obtain a formula for it, see ()5ip below. The \sq) are connected 
to the creation operators according to 

\sq) = ^^Xr,) ■ ■ ■ tp\x,)a\yj ■ ■ ■ a\y,m , (46) 

where |0) G denotes the vacuum state. A relevant fact is that the creation and annihi- 
lation operators ip^ipja^ and a possess kernels. Using the canonical ( ant i-) commutation 



20 



relations for ip and a, one obtains from pUj) the following formulas for the kernels of 
ip{r) and a(r), r G M.^: 

{sqmr)\sq) = 5„,n'-i Sm,m' U r - x') ^^-(^ _ (47) 

(sg|a(r)|sg') = - x') (-l)^(^''^'') 5'"^' (y U r - y') (48) 

where (x, y) = q = vr(g), and x') denotes the permutation that carries x to x' given 
that X = x'. The corresponding formulas for ip"^ and can be obtained by exchanging 
q and q' on the right hand sides of ^7\i and ()48|) . For the smeared-out operator a^(r'), 
we obtain 



We make use of the resolution of the identity 



I = J dq\sq){sq\ . 
Q 

Inserting (jSUj) twice into (jl^j) and exploiting (jTTj) and (HH)), we find 

{sq\Hj\sq') = 5^-1,^' 5^^{x - x') (-l)^(-'-') ^ ^^rn'^^ \ ^ _ ^ ^(^ _ a;) 



(50) 



By fHHj) . the jump rates are 

I Im (^f I sg) (sg | if / 1 sq') {sq' | 



{^\sq'){sq'\^) 
More explicitly, we obtain from (j3T|) the rates 



(51) 



(52) 



(T{q\q) = 6nn' Sm-l,m' ^^"(x - x') ^1 (?/ \ ?/ " ^crea(g' U y|g') 

yes/ 

+ 5nn' 5'"(X - X') J2 ^'"'iy - V \ y') ^ann(g' \ y'W) 



(53) 



with 



Crcrea(g' U 2/|g') 



C^ann(g \ ?/1g') 



2v/m' + 1 



y'ey' 



Im^*(g) (-l)^(^'^') ^{y - x')-^{q') 



h 



^*(g') ^(g') 



hym' 



Im^*(g) (-1)^^(*'^') Yl Vi.y' -x')'^{q') 

x'dx' 



(54a) 
(54b) 
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for arbitrary g' G vr~^(g') and q G vr~^(g) with q = {x', y'Uy) respectively q = (x', y'\y')- 
(Note that a sum sign can be drawn out of the plus function if the terms have disjoint 
supports.) 

Equation ()53p is worth looking at closely: One can read off that the only possible 
jumps are {x',y') {x',y' U y), creation of a boson, and {x',y') {x',y' \ y'), an- 
nihilation of a boson. In particular, while one particle is created or annihilated, the 
other particles do not move. The process that we considered in |[14j consists of pieces of 
Bohmian trajectories interrupted by jumps with rates (|^: the process is thus an ex- 
ample of the jump rate formula (P^ . and an example of combining jumps and Bohmian 
motion by means of process additivity. 

The example shows how, for other QFTs, the jump rates (j^Uj) can be applied to 
relevant interaction Hamiltonians: If Hj is, in the position representation, a polynomial 
in the creation and annihilation operators, then it possesses a kernel on the relevant 
configuration space. A cut-off (implemented here by smearing out the creation and 
annihilation operators) needs to be introduced to make Hj a well-defined operator on 
L\ 

If, in some QFT, the particle number operator is not conserved, jumps between the 
sectors of configuration space are inevitable for an equivariant process. And, indeed, 
when Hj does not commute with the particle number operator (as is usually the case), 
jumps can occur that change the number of particles. Often, Hj contains only off- 
diagonal terms with respect to the particle number; then every jump will change the 
particle number. This is precisely what happens in the model of [Tij . 

3.2 Efficient Calculation of Rates in the Previous Example 

We would like to give another, refined way of calculating the explicit jump rates ()53|) 
from the definition of Hj. The calculation above is rather cumbersome, partly 
because of all the 5's. It is also striking that only very few transitions q' ^ q are 
actually possible, which suggests that it is unnecessary to write down a formula for the 
kernel {q\Hi\q') valid for all pairs g, q' . Rather than writing down all the 5 terms as in 
fl^ . it is easier to specify the possible transitions q' ^ q and to write down the rates, 
such as ()54a|) and ()54b|) . only for these transitions. Thus, for a more efficient calculation 
of the rates, it is advisable to first determine the possible transitions, and then we need 
keep track only of the corresponding kernel elements. 

3.2.1 A Diagram Notation 

To formulate this more efficient strategy, it is helpful to regard \E' as a cross-section of a 
fiber bundle E over the Riemannian manifold Q, or of a countable union E = {j^ E^^ of 
bundles i?*^*-* over Riemannian manifolds Q*^*-* with Q = {_]^ Q*^*^ (In the present example, 
with Q given by (j4U|) . we take i to be the pair (n,m) of particle numbers, Q^"'™-) to be 
the (n, m)-particle sector, and E^^^ to be defined by (|HS|) (with vr the natural projection 
from Q \ A, with Q given by (|^T|) . to Q). The q G vr~^(g) can be viewed as defining an 
orthonormal basis of Eg.) 
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A key element of the strategy is a special diagram notation for operators. The 
operators we have in mind are Hj and its building blocks, the field operators. The 
strategy will start with the diagrams for the field operators, and obtain from them a 
diagram for Hj. The diagram will specify, for an operator O, what the kernel of O is, 
while leaving out parts of the kernel that are zero. So let us assume that O has kernel 
{q\0\q'), i.e., (0^)(g) = ^ {q\0\q') ^ {q') dq' . The diagram 

g' ^ F{q\ A) (55) 

means that the operator O has kernel constructed from F and K, 

{q\0\q') = [ dX6{q-F{q\X))K{q',X), (56) 



where A varies in some parameter space A, F : Q x A ^ Q, and K is a function (or 
distribution) of q' and A such that K{q', A) : Eg' Epi^qi^x) is a C-linear mapping. 

The role of A is to parametrize the possible transitions; e.g., for the boson creation 
()54ap in the previous section, A would be the position y of the new boson, and A = M^. 
The notation does not explicitly mention what A and the measure d\ are; this will 
usually be clear from the context of the diagram. The measure d\ will usually be a 
uniform distribution over the parameter space A, such as Lebesgue measure if A = M*^ 
or the counting measure if A is finite or countably infinite. We may also allow having a 
different Aqi for every q' . 

In words, may be read as: "According to O, the possible transitions from q' are 
to F{q' , A), and are associated with the amplitudes K{q' , A)." In fact, when O = H, a 
jump from q' can lead only to those g's for which q = F{q', A) for some value of A, and 
the corresponding jump rate is 

V V>i, Jl'dj ^*(g')^(g') ' ^ ^ 

provided that for given q', F{q', ■ ) is an injective mapping. Here, cr(g|g') is the density 
of the measure o"((ig|g') with respect to the measure on Q 



Hidq) = j dX6{q-F{q',X))dq, 



(5J 



where S{q — qo) dq denotes the measure on Q with total weight 1 concentrated at go- (ISHl), 
the image of d\ under the map F{q', ■ ), is concentrated on the set {F{q', A) : A G A} 
of possible destinations and plays the role of the "uniform distribution" over this set. 
In other words, ()57|) is the rate of occurrence, with respect to dX, of the transition 
corresponding to A. (For the boson creation rate ()54a|) . fiq'{dq) turns out the Lebesgue 
measure in y on the subset {q' U y : y E M.^ \ q'} Q.) 

Given O, the choice of A, F, and K is not unique. One could always choose A = Q, 
F{q',q) = q, and K{q',q) = {q\0\q'), which of course would mean to miss the point 
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of this notation. The case that F and K do not depend on a parameter A is formally 
contained in the scheme by taking A to be a one-point set (and dX the counting 
measure); in this case means 

{q\0\q') = 6{q-F{q'))K{q'). (59) 

Conversely, whenever = 1, the dependence of F and K on the parameter A is 
irrelevant. 

A basic advantage of the notation (|55|). compared to writing down a formula for 
{q\0\q'), is that many 6 factors become unnecessary. For example, if O is multiplication 
by V{q), then (A is a one-point set and) we have the diagram 

q -^q. 



3.2.2 Operations With Diagrams 

For the product O2O1 of two operators given by diagrams, we have the diagram 

q — ^ F2(Fi(g, Ai), A2) (60) 

O2O1 

with parameter space Ai x A2, for which we also write 

q — ^ Fi[q , Ai) > F2[Fi[q , Ai), A2). (61) 

Ci U2 

We thus define the concatenation of two diagrams by means of the composition of the 
transition mappings and the product of the amplitudes, i.e., using obvious notation, 

qi ^ q2^ qz means qi — > qs. (62) 

Thus, multiplication of operators corresponds to concatenation of diagrams. 

For the sum Oi + O2 of two operators given by diagrams with the same parameter 
space Ai = A2 = A and the same transition mapping Fi{q',X) = F2{q',\) = F{q',X), 
we have the diagram 

, K,{q',X)+K2{q',X) j;.. , .X 

q — — >F{q,X). (63) 

U1+U2 

3.2.3 Diagrams of Creation and Annihilation Operators 

We now write down diagrams for creation and annihilation operators. In the case that 
O = 0{r) arises from formally evaluating an operator- valued distribution 0{x) at x = r, 
the dependence of K{q', A) on A is in the sense of distributions rather than functions. 
More precisely, we have 

K{q',X) = D{q',X)Ko{q\X) (64) 

where D is a (real-valued) distribution on Q x A, and Kq a mapping-valued function 
such that for every q' and A, K^i^q', A) is a linear mapping Egi — > Ep(^q>^x)- 
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For ?/'''" (r) and ip{r), r G M^, we have (recall that x' is a finite subset of 



p3\ 



(x',y')^ix'Ur,y') (#A = 1) (65a) 

{x\y') {x'\x\y') {X = x',\ = x') (65b) 

using linear mappings ctf : Eqi Ei^^'ur^y') ( "append a fermion" ) andef : Eqi E(^x'\x',y') 
("erase a fermion"), which can be regarded as the natural mappings between these fiber 
spaces. They are defined through the following properties: 

af\l/ is appropriately symmetrized (66a) 

(«f^)((x',r),y') = ^L^^{x',y') (66b) 
Vn' + 1 

(ef^) (x, y') = vV^((x, x'),y) (66c) 

where \1/ G £"5', and x is an arbitrary ordering of the set x = x' \ x'. (Recall that the 
set 7r~^(g') of the possible orderings of q' forms a basis of Eg/, so that every ordering 
{x',y') = q' E 7r^^(g') corresponds to a particular component of Thus, {{x',r),y') G 
7i~^{x' U r, y') corresponds to a particular component in E(_x'ur,y')-) 

For the smeared-out creation and annihilation operators a}p{r) and a<^(r'), we have 

(3;',^')Z(^Z!l^(^',^'Uy) (A = M^A = y) (67a) 

al,(r) 

(^j^'^y'^ I^y_^:lf^ (^jr'^y'\y') [\ = y\\ = y') (67b) 

where ctb ( "append a boson" ) and ( "erase a boson" ) are the analogous linear mappings 
relating different spaces, ctb : Eqi — > -E(x',y'uy) and ^b : -Eg' — * ^{x',y'\y'), defined by the 
following properties: 

ab^ is appropriately symmetrized (68a) 
(«b*)(x',(y',?/)) = ^2==m{x\y') (68b) 
(£b^')(a;',y) = ^/^^>{x',{y,y')), (68c) 
where y is an arbitrary ordering of the set y = y' \y', x' one of x' , y' one of y' , and 

^ G Eg/. 

3.2.4 Application of the Diagram Method 

Now let us apply the strategy to the example of the previous section. For 

%p'^{r) a]^{r) ip{r), we have the diagram 

q > [x \x ,y) > [x \ x ,y U y) > {x \x Ur,y U y) 

V-W at(r) i'Hr) 
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with A = x' X M'^. Using the concatenation rule ()62|) . we can write instead 



{x' \x'Ur,y'U y). 



Integrating over dr, we obtain, since x' \ a;' U r may be replaced by x', which is inde- 
pendent of x\ 



'P{y-x')aia^ei 

q > [x ,y Uy), 



, ^ (69) 

J dr^p'l{r) alp{r)tli(r) 

with A = R'^. We have now taken care of one of two terms in involving rather 
than a. From (jH^ll we read off, without a big calculation, that this term corresponds 
to jumps (x', y') (x', y' Uy), or creation of a boson. The corresponding jump rate is 
given by (jSZj), and reads here: 



a{x\y'Uy\q') = - 



Im^*(x',?/'Uy) Yl Viy-x')aia^ei'^{q') 

x' dx' 



(70) 



This result agrees with ()54a|) .^ 

We treat the term J dr (r) a^ir) ipif) in the same way: We begin with the diagram 

q > (x \x ,y) -— ^ (x \x ,y \ y) — — > (x \x Ur,y \ y) 

■0(r) a¥>{^) ipi{r) 

with A = x' X y'. Then we integrate over dr and obtain the associated jump rate 



a{x,y \y\q) 



lm<i/*{x',y'\y') ^ ip{y' - x') afe^Ei^iq') 
x'ex' 



(71) 



h ^*(g')^(g') 

which agrees with ()54b|) . Finally, Hj (the sum of both contributions) corresponds ac- 
cording to (f^^ to jumps which, since the two contributions have no transitions q' —>■ q 
in common (or, in other words, since their kernels have disjoint supports in Q x Q), are 
either q' — > {x',y' U y), with rate (fTOIl . or q' — > {x',y' \ y'), with rate (f7T|) . 



3.3 Pair Creation in an External Field 

As our second example, we present the Bell-type version of a reasonable and often used 
QFT of electrons and positrons, in which the electromagnetic field is a background field 
|31j . The Bell-type version exhibits pair creation and annihilation (in the literal sense) 
and employs various notions we have introduced: process additivity, the configuration 
space r^(M^) of a variable number of identical particles, the free process, POVMs which 
are not PVMs, and stochastic jumps. 

^Here is why: First, '^*{q')'i'{q') — n'! m'! ^'(g') because the inner product in Eq/ involves 

summation over all q' G Tr^^{q'). Similarly, the square bracket in the numerator of H7U|) involves the 
inner product of -E(a:'.j^'uy') ' consisting of n'l (m' + 1)! contributions. The numberings q and q' in H54a|l 
can be so chosen that x = x' , x' gets the last place of x' , and y = y' U y'; then g{x, x') is trivial, and 
afabEf^l?) = (n')~^^^(?7?/ + iy^/'^{n'Y/'^'^{q'). Thus, the square bracket in is n'l m'\y/m' + 1 
times the square bracket in (|54a|) . 
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3.3.1 Fock Space and Hamiltonian 

We consider the second quantized Dirac field in an electromagnetic background field 
Afj^{x,t). In terms of field operators, the Hamiltonian reads 

H= J d^x : <i>*{x)[-icha-V + pmc^ + e{cf A + Ao)]<i>{x) : , (72) 

with colons denoting normal ordering. Note that H is time-dependent due to the time- 
dependence of A^{x,t); more precisely, Hj is time-dependent while Hq is fixed. As a 
consequence, the relevant jump rate (j29p is now time- dependent in three ways: through 
Hj, through \E', and through q' = Qt- 

We quickly recall what the Hilbert space and the field operators are, and specify 
what POVM we use. After that, we construct the associated process. 

The Hilbert space L^(M^,C^) of the Dirac equation is split into the orthogonal sum 
© of the positive and negative energy subspaces of the free Dirac operator, 

= —ichcx ■ V -|- /3mc^ . 

The 1-electron Hilbert space and the 1-positron Hilbert space are copies of 
and the Fock space ^ = TJi^^^^ arises then from the one-particle Hilbert space 
^(1) = ^ © ^ in the usual manner: with the anti-symmetrization operator Anti , 

CXD 

^ = 0Anti((jre©^p)®^), (73) 

N=0 

which can be naturally identified with 

oo oo 

:= ® = Anti (^e®") ® Anti (^p®^) . (74) 

n=0 n=0 

Since C L^(]R^, C''), can be understood as a subspace of 

oo oo 

=^ext := Anti (L2(M^ C^)®") ®0 Anti (L2(M^C^)«"). (75) 

n=0 n=0 

We choose the POVM and configuration space in the way suggested by the form 
dZH), rather than 

Q = T^iR') X T^R'), (76) 

where the first factor represents electrons and the second positrons. (Recall from Sec- 
tion Em that T^(R^) denotes the space of all finite subsets of R^. Another interesting 
possibility, suggested by the representation (f73|) . is to set Q = r^(M^). This would 
mean that, insofar as the configuration is concerned, electrons and positrons are not 
distinguished. However, we will not pursue this possibility here.) The natural POVM P 



27 



(see Section 14.2.31 and Section I2.11|) can be expressed as an extension from rectangular 
sets (the existence of such an extension is proved in Section 4.4 of [15 J: 

p{B, X Bp) = rp(^)(5e) ® rp(i)(5p) 

with P(^) the POVM on that we considered before, arising by projection from the 
natural PVM on L^(R^,C'^). Alternatively, P can be viewed as arising, by projection 
to je, and from Q = U^r=o(^^)'' ^ U^=o(^^)^ ^o Q, of the natural PVM on Q acting 
on Jifext- Note that P represents the usual |\E'p distribution in the sense that for a 
configuration q with electrons at Xi, . . . ,Xn and positrons at 5i, . . . , 5^, we have 

P(rfg) = (^lP(rfg)l^) = n\n\\¥''^^\xi,...,xn)\'^dxi---dxn 

where is just the wave function (]R^)"+" — ^ (^(Q4^(g)(n+n) when we decompose 

the state vector in the manner suggested by (|75|) . \E' is normalized so that 

oo „ 

I dx, - ■ ■dXn\¥''~^\xu . . . ,XJ,)\^ = 1. 

n,n=0 

The field operator is defined by 

$(/)=6(P+/)+rf*(CP_/) (77) 

where / is a test function from ^^(^3^ C^), P± is the projection to C L^{R^, C^), C 
is the charge conjugation operator which maps to Jif+ and vice versa, and b is the 
electron annihilation and d* the positron creation operator. Letting be the standard 
orthonormal basis of C^, i = 1,2, 3, 4, $(a;) stands for ^i{x) = $(ej 6{ ■ — x)), where i 
gets contracted with the a matrices. Similarly, we define, as usual, 

h{x)=b(^P+{e,6{- -x))) (78a) 
and di{x) =d(CP-{e,6{- -x))y (78b) 

We thus have ^i{x) = bi{x) + d*{x). 



3.3.2 The Associated Process 

We now describe the associated Markov process. The free part of (|7^ . 

preserves particle numbers (it commutes with the electron and positron number opera- 
tors), evolving the (n, n)-particle sector of the Fock space according to the free (n, n)- 
particle Hamiltonian 

k=i 
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with 



where o:*^'^-' and /^'-^^ act on the fc-th electron index in the tensor product representation 
(|74p and S*-'^-' and P^'^^ on the k-th positron index. is the gradient with respect to 

With Hq is associated a deterministic motion of the configuration in Q, the free pro- 
cess introduced in Section ITHl During this motion, the actual numbers^ A^, N of^electrons 
and positrons remain constant, while the positions (Xi, . . . , X^, JCi, . . . , X^) =: Q 
move according to Bohm-Dirac velocities (j31|) . i.e. 

. vi/-(g)aWvi/(g) 



^*(Q)^(Q) 
^ *(g)5(~')^(Q ) 

^*(Q)^(Q) 



where numerators and denominators are scalar products in (C^)*^*^^"*"^-*. 

We turn now to the interaction part. Setting A = a ■ eA + cAq, we have that 

Hi = J d^x : $*(a;) A{x) <l>{x) : = (80a) 
= J2 I d'^x : {b*{x) + di{x)) A'^'ix) {b,{x) + d*{x)) : = (80b) 

4 ^ 

= V / rf=^a; f6*(a;) A''^{x) bj{x) + d,(a;) A*'^(a;) bj{x) + 

..=1 ^ (80c) 

+ b*{x) A''^{x) d*{x) - d*{x) A''^{x) d^{x)y 

Since Hj is a polynomial in creation and annihilation operators, it possesses a kernel and 
corresponds to stochastic jumps. To compute the rates, we apply the strategy developed 
in Section 13.21 using diagrams. To this end, we regard fermionic wave functions again 
as cross-sections of a bundle E, defined here by 

= {C^f^^iCY^. (81) 
Fermionic symmetry of a cross-section of means that 

£>(il...in),5(*l---*ri) il . . .i„ ,11 . . .Ifi 

^{g{Xi...Xn),d{Xl...Xn)) = {-If (-1)^ ^ {Xi . . . Xn,Xi . . . X^) (82) 
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for all permutations g & Sn and g & Sn. 

The diagrams for b*{x),bi{x), d*{x), and di{x) are 



(x', x') > (x' U x', x ) (83a) 









S+^(a3'-x) £e(ej) 




> 

bi{x) 








d*{x) 




S-l(x'-x)e^(e,) 


di(3;) 



^{x'\x\x') (83b) 
, X U X ) (83c) 



di(3;) 

where the matrix function 5*+* (a:;) is defined as the j-component of P+(ej(5(-)), and 
S-j{x) as the j-component of CP_(ej5(-)). The linear mappings ae(ej) : -Eg' — > 
-E{a;'ua;',x') ("append an electron with spinor ey) and Se^Bj) : E(^x'\x',x') ("erase 

an electron, contracting with spinor e^") are defined through their properties that for 

a^^ is appropriately symmetrized (84a) 

iae{ej)m){{x,x'),x') = , ^ ^{x^,x')®ej (84b) 

\/ri! + 1 

(£e(ej)^)(x,5^') = v^^j((x,cc'),x'), (84c) 

where x is an arbitrary ordering oi x = x' \ x' ^ x! one of x', and x' one of x'. We refer 
to the last electron slot when writing the tensor product or taking the j-component. 
ap(ej) and ep(ej) are defined analogously. 

For the four terms in ()80cp . we thus get the four diagrams (omitting the multiplication 

by A^^\x)) 

hjix) b*{x) 
, ^, EkS+']ix'-x)e4^,) , E,5_f(i'-x)ep(e,) . 

(x , X j > [x \ x ,x ) > [x \x ,x \ x ) (85bj 

bj{x) di{x) 
, ^, EkS^']ix'-x)a^ie,) , ^, , Z,S+Ux'-x)a4e,) ^ ^ ^ , 

(x , X ) > (x , X U a:; ) > (x Ux ,x Ux) (85c) 

d*{x) b*{x) 

(x , X j > (x , X \ x ) > (x , X \ x U a; ). (85dj 

We read off that the first term corresponds to the jump of a single electron from x' to 
x" , while all other particles remain where they were, the second to the annihilation of 
an electron-positron pair at locations x' and x\ the third to the creation of an electron- 
positron pair at locations x' and 5', and the last to the jump of a positron from x to 
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x' . The corresponding jump rates are 

ae(x \ x\^x ,x|g) = ^(^.^ (86a) 

, ^ ^ [(2M)Imv^*(g)EMxL'n(^',2')^p(e.K(e.)vl/(g')]+ _ , 
aa„n(x \ a; ,x \ 03 |g ) = ^.(^.^^(g,) (86b) 

^ [(2/n)ImvI/*(g)^,,Xc'rL(^',5')«e(e,)ap(e,)vI/(g')]+ , 
a_(a: Ua.,x Ua.|g) = ^^MW^ ^ ^ 

a^{x,x\xV^x\q) = ^^*{(^)^{(^) ' ^^^^^ 

where q denotes the respective destination, and 

X,\i(a;', 5') = E / ^-(^' - ^) ^'''(^) ^4(^' - ^) (87b) 

X'.L{^\ 5') = E / ^+'(^' - ^) ^''(^) ^-'(^' - ^) (87c) 
xY{x\ x") = -J2 fd'x S.]{x" - x) A^'^ix) S4{x' - x). (87d) 

The process for Hq + iJ/ that we obtain through process additivity is the motion (f?^ 
interrupted by stochastic jumps with rates (jHUj) . 

Note that the jump of a single electron has small probability to be across a distance 
much larger than the width of the functions S±, which is of the order of the Compton 
wavelength of the electron. Similarly, the distance |cc — 5| of a newly created pair, or 
of a pair at the moment of annihilation, has small probability to be much larger than 
the width of 5"^. While the jump of a single electron or positron leaves the number N 
of electrons and the number of positrons unchanged, pair creation and annihilation 
can only either decrease or increase both and by 1. As a consequence, the actual 
net charge N — N is conserved by the process. 



4 Second Quantization of a Markov Process 

4.1 Preliminaries Concerning the Conditional Density Matrix 

In the next section, we describe the algorithm for the "second quantization" of a process. 
But before that, we have to introduce, as a preparation, the notion of a conditional 
density matrix. In jTHj, we have defined for Bohmian mechanics the conditional wave 
function of, say, subsystem 1 of a composite system with configuration space Q = 
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Qi X Q2 by ^cond(Q'i) = ^2)- From a complex wave function \1/ : Q C, together 

with the actual configuration Q2 of the environment of the subsystem in the composite, 
we thus form a wave function ^cond : Qi — C; for Bohmian mechanics with spin, in 
contrast, we would not, in general, obtain a suitable wave function for subsystems in 
this way, because ^'cond as just defined would have more spin indices than appropriate. 
We can however still define the conditional density matrix for subsystem 1, 

W,,r.ds^,s'Mi,q'i) = - J]^si,.2(9i,Q2)^':i,,,(gi,Q2) (88) 

where the s's are spin indices. In order that W, like any density matrix, have trace 1, 
the normalizing factor 7 must be chosen as 



This W can play most of the roles of the conditional wave function in spinless Bohmian 
mechanics. The notion of a conditional density matrix easily generalizes from the sit- 
uation just described, corresponding to wave functions in L^(Q,C'^) and the natural 
localization PVM, to the situation of any product localization POVM on any tensor 
product Hilbert space: for = M{ ® ^ and P{dq\ x dq2) = Pi{dqi) (8) P2{dq2), set 



tr2{\^){^\ PiQiXdq2)) 

l^cond ~ 



tT{\^){^\PiQlXdq2)) 



q2=Q2 



(89) 



where tr2 is the partial trace over The quotient is to be understood as a Radon- 
Nikodym derivative in q2. Like conditional wave functions, conditional density matrices 
cannot be defined in orthodox quantum theory, for lack of the configuration Q2. We 
stress that conditional density matrices have nothing, absolutely nothing, to do with 
statistical ensembles of state vectors in J^. Like any density matrix, they do, however, 
define a probability distribution on Qi, 

pfcond(.)^tr(iyeondPl(-)), (90) 

which coincides with the conditional distribution of Qi given Q2, 

[^\Pi{-)(»P2{dq2)\^) 



P(<5i e • IQ2) = 



P2{dq2)\^) 



The evolution of PFcond is not autonomous; it will typically depend on (and always 
be determined by) and Q2,t- For a given density matrix 14^ of a system that is 
not regarded as a subsystem, however, one can define (as usual) the time evolution 
by Wt = e~*^*/^ W e*^*/'^, which gives rise to a time-dependent distribution P'^* ( " ) ~ 
tr( WtP( • )). We call a Markov process that is P^*-distributed at every time t equivariant 
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with respect to W and H. Given the right initial distribution, this is equivalent to the 
following condition on the generator: 



='^lmtT{WP{-)H) . (91) 

This is the version of ()12j) for density matrices, and defines an equivariant generator 
with respect to W and H. 

Since conditional density matrices will play a crucial role in the construction of the 
many-particle process, we require that, as part of the input data of the algorithm, we 
are given an equivariant generator =^14^^ for every density matrix from a dense subset 
of the density matrices in Jif^^^* (g) J^(^) . This is not much of a restriction, as all rele- 
vant examples of equivariant generators naturally extend to density matrices: Bohmian 
mechanics with spin space C*"' can be extended jl| to 

\ fcT ^gt^c* ^(?, / / N fao\ 



Bohm-Dirac to 



and minimal jump rates to 



^ _ trc4(iy(g,g)a) 



Note also that (jH^ would not make any sense if W represented a statistical ensemble |3], 
whereas it makes good sense for conditional density matrices, expressing the true relation 
between the Bohmian velocity for a subsystem arising from and the conditional 
density matrix (|HHj) of that subsystem. Mutatis mutandis, the same is true of (|^. 
Similarly, in case that P is a PVM, expresses the jump rates for a decoupled 
subsystem arising from (j29|) for the composite in terms of the conditional density matrix 
of that subsystem. 



4.2 Algorithm 

The input data of this algorithm are the one-particle Hilbert space J^^^\ configuration 
space POVM P^^\ and a family of generators .Sf-^^ = labeled by the density 
matrices W from a dense subset of the density matrices in J^^^^* (g) J^f^^^ . The output 
is a family of generators TJif^^^ = =Sfo = -^o,* labeled by the state vectors \l/ in (a dense 
subspace of) Fock space. If is equivariant with respect to W and then =Sfo,* 
is equivariant with respect to and Hq. 

The algorithm is based on two procedures for suitably combining generators for direct 
sums or tensor products of Hilbert spaces. 
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4.2.1 Direct Sums 

Given a finite or countable sequence of Hilbert spaces witli POVMs P(") on con- 

figuration spaces Q^"^ , and for eacli n a family of generators labeled by the vectors 
in J^^^\ there is a canonically constructed family of generators =Sf® = J^^, labeled 
by the vectors in the direct sum The space Q in which the corresponding 

process takes place is the disjoint union of the Q^'^l If every ^^^J is equivariant with 
respect to \E'„ G Jif'^'^^ and H'^'^\ then is equivariant with respect to \E' G 0„^*'"'' 
and ©„i/("). 

Here are the details. The POVM P = 0„ p(") on Q that naturally arises from the 
data is given by P{B) = ®^P^''\B n Q^")) for C Q. Let Pn denote the projection 
^ — >• The generator ^® is given by 

P)|g{n) = =^P„i/||P„-I.|j (p|Q{n)) • (95) 

It generates a (Markov) process Qf such that when Q'^ G Q*-"-*, it is generated by 
the state vector P„\E'/||Pji\E'||, i.e., it is a Markov process qI"'' in Q^"^ generated by 
-^p"^/IIP^^I|- The equivariance statement follows directly, since ||P„\E'j||^ = Pt(Q^"'') is 
invariant under the evolution generated by Hq = ©^-ff^"-*- 



4.2.2 Tensor Products 

Given a finite sequence of Hilbert spaces J^'^l, . . . , ^'"1 with POVMs P^ on con- 
figuration spaces QW, and for each i a family of generators .Sf'*' = labeled by 
the density matrices on there is a canonically constructed family of generators 

<if® = labeled by the density matrices on the tensor product ^'^1 ® ■ ■ ■ ® ^'"l 
The corresponding process takes place in the Cartesian product Q = Q'^^ x ■ ■ ■ x Qt""'. 
If every is equivariant with respect to the density matrix Wi on the 
Hamiltonian ifW, then is equivariant with respect to W on ® ■ • ■ ® ^t""' and 
= E 1 ® ■ • • ® ifH ® ■ • • ® 1 = ^ i/i. 

i i 

Here are the details. The POVM that naturally arises from the data is^° 

P{dq, X • ■ • X dqj = P^^\dq,) ® • ■ ■ ® P^^'^dqJ. (96) 

For any q E Q, let denote its i-th component and let = (qi, . . . , qi_i, qi+i, ■ ■ ■ , qn)- 
For every i and g^, define 

^ ti^,{WP{dq,x...xQl^x...xdqJ) 
'^^'^ tT{WP{dq^ X • ■ ■ X X ■ ■ ■ X rfq„)) ' 

where tr^j is the partial trace over all factors except J^^W. This Wi is the conditional 
density matrix, regarded as a function of the configuration % of the other particles. Now 

^°The existence of the tensor product POVM is a consequence of Corollary 7 in Section 4.4 of |15| . 



34 



consider the process on Q according to which the i-th particle moves as prescribed by 
<5f^^. while the other particles remain fixed. The generator of this process is 

pMdqi) (97) 

where p^i is the marginal distribution of Qi (i.e., p integrated over q^) and p{ ■ |gj) is the 
conditional distribution of given Qi = the square bracket is a function of qi and 
a measure in dq^. Now define .^^p = '^^ip. 

i 

To see that =2"^ is equivariant when the =Sft*' are, we have to check ()9H) . Note first 
that P^(rfqjgi) = iT{Wi{qi) P^^{dq^)). Due to the equivariance of ^W, for p = the 
square bracket in (j97p equals (2/^) Imtr(VFj(gj) pW((iqJ if W) , from which we obtain 
(1^ for and Hi and hence for ^® and H. 

The definition of =Sf® reproduces the many-particles Bohm law (0) with or without 
spin from the one-particle version (or, for distinguishable particles, from several different 
one-particle versions having different masses and spins). Similarly, it reproduces the 
many-particles Bohm-Dirac law (j!?T|l from the one-particle version. 

4.2.3 Second Quantization of the POVM 

Let denote the space of all subsets-with-multiplicities of Q^^^ having n elements 
(counting in the multiplicities). P^^^ naturally defines a POVM p(^)®"' on (Q'-^-')" acting 
on ^(1)®'^ by P(i)®"(c/qi x ■ ■ ■ x dq^) = P^^\dq^) ® ■ ■ ■ ® P^^\dq^), and a POVM 
p("-) on Q*^") acting on = Pj_ (the n-particle sector of Fock space, with P± 

the projection to the subspace of (anti-)symmetric elements of ^(^)®"^ depending on 
whether we deal with fermions or bosons) by 

p(n)(5) = P(^)«"{(gi, . . . , e {Q^'^r ■■ ...,qn}eB} 

for B C Q^"\ where {qr^,...,q„} should be understood as a set-with-multiplicities.^^ 
Since P^^^B) is invariant under permutations, it maps symmetric to symmetric and anti- 
symmetric to anti-symmetric elements of jifW'^" and thus acts on ^("^ for bosonic or 
fermionic Fock space. The corresponding POVM on Q is then P = FP'-^^ = ©^P*^"^; 
more precisely, for B C 

oo 

p(p) = 0pW(pn q(")). 

n=0 

^^This agrees with the definition given in Section [3.11 for the case of a PVM and the coincidence 
configurations removed from configuration space. 
i^In case that P'-^^ is nonatomic, 

pin) 

can equivalently be defined in the following way: For the set A 
of coincidence configurations we set P*^"^ (A) — 0, and for volumes dq^, . . . , rf<7„ in Q'^^^ that are pairwise 
disjoint, we have a corresponding volume dq in Q^"\ which can be obtained from dqi x • • • x (ig„ C 
(g(i))« by forgetting the ordering, and we set P^"'){dq) ^ nl P± P^^^\dq^) (g) ■ ■ ■ (g> P^^\dq^) P±. 



^^p■.-- 



35 



4.2.4 Construction of the Free Process 

Equipped with the two procedures for direct sums and tensor products, we complete 
the construction of the free process. 

The "tensor product" procedure above provides a process on (Q'-^-')" from n identical 
copies of J!f^^\ For a state vector ^f*^"^ e ^(") = P^^(^)®" from either the symmetric 
or the anti-symmetric elements of the n-fold tensor product space, let W be the pro- 
jection to \E'("-); the generator is permutation invariant because the tensor-product 
construction of is permutation covariant and a permutation can at most change the 
state vector by a minus sign, which does not affect the density matrix. Consequently, 
the ordering of the configuration is irrelevant and may be ignored. We thus obtain a 
process on Q^"'^ whose generator we call We now apply the "direct sum" procedure 

to obtain a process on Q. 

5 Towards a Notion of Minimal Process 

In this section, we investigate the common traits of the Markov processes relevant to 
Bell-type QFT, which can be summarized in the notion of a minimal process associated 
with "^,11, and P. We begin with a closer study of the minimal free generator (jH^ . and 
then explain why we call the minimal jump rates "minimal." Finally, in Section 1^31 we 
give an outlook on the notion of minimal process. 

5.1 Free Process From Differential Operators 

In this section, we discuss some of the details, concerning the two equivalent formulas 
fl32|) and (jH^I) for the backward and forward version of the minimal free generator in 
terms of H, P, and that we omitted in Section 12.91 To begin with, L as defined 
by (j32p satisfies some necessary conditions for being a backward generator: Lf{q) is 
real, and LI = where 1 is the constant 1 function (this corresponds to ^p{Q) = 0, 
or conservation of total probability). In case L is indeed a backward generator, the 
corresponding process is equivariant because 

^¥{dq) ^ Re{^\il[H,P{dq)]\^) = hm {^\P{dq)H\^) ^ F{dq) . 
h n 

One way to arrive at formula (j32j) has been described in Section 12.91 A dif- 
ferent way, leading to (jHS)), is to start from the ansatz =Sfp = where A de- 
notes a (signed-measure-valued) linear operator acting on functions. Equivariance 
means Al{dq) = {^\{[H, P{dq)]\^) . This suggests Af{dq) = {^\f {[H, P{dq)]\-^), or 
Af(dq) = (\E'| I [if, P((ig)] /|\&), or a convex combination thereof. Since Af{dq) must be 
real, we are forced to choose the combination with coefficients | and |, or equivalently 
Af{dq) = Re{^\f j:[H,P{dq)]\-^), which is 

That ^ generates a deterministic process (when it is a generator at all) is suggested 
by the following consideration — at least when H and P are time-reversal invariant: 
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replacing \l/ in by T\l/ where T is the anti-hnear time reversal operator (see Section 
16. 1|) changes the sign of ^ . The only generators ^ such that — =Sf is also a generator 
are, presumably, those corresponding to deterministic motion. 

This gives us an opportunity to check for which H fj32|) does define a process: for 
a deterministic process we must have L = f ■ V where v is the velocity vector field. 
It is known that vector fields, understood as first-order differential operators, are those 
linear operators L on the space of smooth functions that satisfy the Leibniz rule L{fg) = 
f'Lg + gLf. (j32p is certainly linear in /, so we have to check the Leibniz rule to see 
whether L is indeed of the form v ■ V and thus the backward generator of a process. 

We can see no reason why L should satisfy a Leibniz rule unless P is a PVM, which 
implies that 

fP{dq) = f{q)P{dq), (98) 
and H is such that for all (nice) functions / and g, 

[[HJlg]=h (99) 

for some function h, which holds if is a differential operator of order < 2. (If H = —A, 
then h = — 2V/ ■ V^^; if i7 = —ia V for whatever vector of matrices ex., or if H 
is a multiplication operator, then h = 0.) To check that the Leibniz rule is obeyed 
in this case, note that we then have that [H,fg] = [H,fg] = [H,f]g + f[H,g] = 
f[H, g] + g[H, /] + [[H, f],g] . Using this in (1221), we find that, due to dMl), the first two 
terms give the Leibniz rule, whereas the last term, due to (j^ . does not contribute to 
the real part in ((221) • 

When is an space over Q and P the natural PVM, i.e., when \l/ is a function, 
(I32|) can be written in the form 

(100) 

'^^^ h ^*(g)^(g) ^ ^ 

where / is the multiplication operator corresponding to /. From this, one easily reads 
off the Bohm velocity (P) for the A^-particle Schrodinger operator (jS} with or without 
spin. Similarly, we get the Bohm-Dirac theory when H is the Dirac operator in = 
Anti L^(M^, C^)®^, Q the manifold of subsets of with elements, and P the obvious 
PVM. (ITnn|l also leads to the Bohm-Dirac motion if = L\R^, C^)®^, Q = R^^, and 
P is the natural PVM, but not if is the positive energy subspace because then the 
appropriate POVM P is no longer a PVM. 

To see that the "second quantization" algorithm maps minimal free generators to 
minimal free generators, or, in other words, preserves the relation ()35j) between Hamil- 
tonian and generator, observe first that ()35|1 naturally extends to density matrices, and 
the extension, if a generator, is equivariant. Next check that the "direct sum" and 
"tensor product" procedures of Section ^21 are compatible with ()35|) when P is a PVM. 
Finally, observe that the (anti-)symmetrization operator commutes with the n-particle 
Hamiltonian, with P{B) for every permutation invariant set B C (Q^^^)", and with / 
for every permutation invariant function / : (Q^^-*)" M. 
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5.2 Minimality 

In this section we explain in what sense the minimal jump rates ()29j] — or ()22j] — are 
minimal. In so doing, we will also explain the significance of the quantity J defined in 
()28|) . and clarify the meaning of the steps taken in Sections 12 . 61 and 12 . 71 to arrive at the 
jump rate formulas. 

Given a Markov process Qt on Q, we define the net probability current jt at time t 
between sets B and B' by 

UB, B') = jim [Probjgt e B' , Qt+At e B} - (101) 



-FToh{QteB,Qt+AteB'} 



This is the amount of probability that flows, per unit time, from B' to B minus the 
amount from B to B'. For a pure jump process, we have that 

Jt{B,B')= j at{B\q') pt{dq') - j at{B'\q) pt{dq) , (102) 

g'eB' q&B 

SO that 

UB,B')=j^,,iBxB') (103) 
where jo-,p is the signed measure, on Q x Q, given by the integrand of (fT3j) . 

j,,p(rfg X dq') = a{dq\q') p{dq') - cr{dq'\q) p{dq) . (104) 

For minimal jump rates a, defined by ()29|) or (j^^ (and with the probabilities p given 
by (j24|l . p = ¥), this agrees with (j28|l . as was noted earlier, 

ja,p = J*,H,P , (105) 

where we have made explicit the fact that J is defined in terms of the quantum entities 
\E', H, and P. Note that both J and the net current j are anti-symmetric, J*'' = — J and 
j*"^ = —j, the latter by construction and the former because H is Hermitian. (Here tr 
indicates the action on measures of the transposition (g, q') (g', q) on Q x Q.) The 
property ()105p is stronger than the equivariance of the rates a, ^a^t = dFt/dt: Since, 

by 

{^^p){dq)=j^,,{dqx Q), (106) 

and, by (j2Hl), 

d¥ 

(dq)=I{dqx Q), (107) 

dt 

the equivariance of the jump rates a amounts to the condition that the marginals of 
both sides of (jl05j) agree, 

Ja,p{dq X Q) = S{dq X Q) . (108) 
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In other words, what is special about processes with rates satisfying ()105|) is that not 
only the single-time distribution but also the current is given by a standard quantum 
theoretical expression in terms of H, ^, and P. That is why we call (jlU5j) the standard- 
current property — defining standard- current rates and standard- current processes. 

Though the standard-current property is stronger than equivariance, it alone does 
not determine the jump rates, as already remarked in j^HSD]- This can perhaps be best 
appreciated as follows: Note that (jl(J4j) expresses jo-,p as twice the anti-symmetric part 
of the (nonnegative) measure 

C{dq X dq') = a{dq\q') p{dq') (109) 

on Q X Q whose right marginal C{Q x dq') is absolutely continuous with respect to p. 
Conversely, from any such measure C the jump rates a can be recovered by forming the 
Radon-Nikodym derivative 

/ , I /X '^(dq X dq') 

pidq') • ^'''^ 

Thus, given p, specifying a is equivalent to specifying such a measure C. 
In terms of C, the standard-current property becomes (with p = P) 

2AntiC = J. (Ill) 

Since (recalling that J = JI+ — J" is anti-symmetric) 

J = 2AntiJ+, (112) 

an obvious solution to ()11H) is 

C = JJ+, 

corresponding to the minimal jump rates. However, (|l(J5p fixes only the anti-symmetric 
part of C. The general solution to (jlllj) is of the form 

C = S+ + S (113) 

where S{dq x dq') is symmetric, since any two solutions to (jlll|) have the same anti- 
symmetric part, and S > 0, since S = C A C^^, because J+ A (J"'')*'' = 0. 
In particular, for any standard-current rates, we have that 

C>r, or a{dq\q)> . (114) 

Thus, among all jump rates consistent with the standard-current property, one choice, 
distinguished by equality in ()114p . has the least frequent jumps, or the smallest amount 
of stochasticity: the minimal rates ()29|). 
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5.3 Minimal Processes 



We have considered in this paper minimal jump processes, i.e., jump processes with rates 
()29|). associated with integral operators H. There is a more general notion of minimal 
process, such that there is a minimal process associated with every Hamiltonian from 
a much wider class than that of integral operators; a class presumably containing all 
Hamiltonians relevant to QFT. This will be discussed in detail in a forthcoming work 

Bohmian mechanics is, in this sense, the minimal process associated with the 
Schrodinger Hamiltonian Q. The minimal process associated with an integral oper- 
ator is the jump process with minimal rates. When the minimal free generator (jH^ 
exists, i.e., when (j32j) is a generator, it generates the minimal process associated with 
H . The minimal process associated with the Hamiltonian of a QFT is the one we have 
obtained in this paper by means of process additivity. The concept of minimal process 
directly provides, perhaps always, the process relevant to a Bell-type QFT. 

To begin to convey the notion of the minimal process, we generalize the standard- 
current property (cf. Section from pure jump processes to general Markov processes: 
the net probability current j of a Markov process defines a bilinear form 

3t{f^g) = hm l-E{f(Q,^^,)g(Q,) - fiQMQt+At)) = {g.Uf) - {f^Ug) (115) 

where Lt is its backward generator, and ( , ) on the right hand side means the scalar 
product of L'^{Q,pt). Then the Markov process satisfies the standard- current property 
if pt = ft and (for / and g real) jt{f, g) is equal to 

Uf.9) = l'^^{'^t\fHg\^t), (116) 

or, in other words, if twice the anti-symmetric part of its backward generator Lt agrees 
with the operator corresponding to Jj as given by (Jt/, g) = 3t{f, g) , 2 Anti Lt = The 
minimal process is then the standard-current process that has, in a suitable sense, the 
smallest amount of randomness. 

Let us consider some examples. The diffusion process with generator ^ given below 
(and for p = P) has the standard-current property (in fact, because its "current velocity" 
is v) for the Schrodinger Hamiltonian Q but is not minimal: 

= ^Ap - div (pv), with v := v + ^V(log (117) 

where A is any positive constant (the diffusion constant) and v is the Bohmian veloc- 
ity (IH); this process was already considered in PUCD]- Note that Nelson's stochastic 
mechanics [21] corresponds to X = h. It is obvious without any mathematical analysis 
that the smallest amount of stochasticity corresponds to absence of diffusion, A = 0, 
which yields Bohmian mechanics. Processes like the diffusion (I117jl for A > seem less 
natural for the fundamental evolution law of a physical theory since they involve greater 
mathematical complexity than is needed for a straightforward association of a process 
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with H and Examples of processes that do not have the standard-current property, 
for the Schrodinger Hamiltonian (j2)), are provided by the alternative velocity formulas 
considered by Deotto and Ghirardi one can say that their current is not the one 
suggested by H and 

We return to the general discussion of the minimal process. As we have already 
indicated, when, for a standard-current process, we view J as well as its backward 
generator L as operators on P), then ijl is the anti-symmetric (skew-adjoint) part 

of L; thus, only the symmetric (self-adjoint) part of L remains at our disposal. Since one 
of the properties of a backward generator is LI = 0, the first possibility L for L that may 
satisfy the formal criteria for being a backward generator is Lf = ^J/ — (^Jl)/. When 
P is a PVM, this is also the operator we obtain by applying, to an arbitrary quantum 
Hamiltonian H, the formula ()32|1 for what we called the minimal free generator, which 
we repeat here for convenience: 

{^\P{dq)\^) ■ ^^^^^ 

Whereas this formula merely provided an alternative definition of the free process in 
Section 12.91 it now plays a different role: a step towards obtaining the minimal process 
from the Hamiltonian H. As we have pointed out in Section 12.91 L is also an obvious 
naive guess for the backward generator L, quite independent of equivariance or the 
current J, since j^[H,f] is the time derivative of /. Moreover, it manifestly satisfies 
LI = 0. For the backward generator L of a standard-current process we must have, 
when P is a PVM, that L = L + S where is a symmetric operator and 5*1 = 0. 
For the minimal process, we have to choose S as small as possible — while keeping 5* 
symmetric and L a backward generator. 

Suppose P is a PVM. Observe then that if is a differential operator (as Hq often is) 
of the kind considered in Section 1231 L is itself a backward generator, so that 5* = is a 
possible, and in fact the smallest, choice. If H is an integral operator, what keeps L, an 
integral operator as well, from being a backward generator is that the off-diagonal part 
of its P-kernel {q,Lq') = F{q)L{q,q') = \lm {^\q){q\H\q'){q'\^) may assume negative 
values whereas the off-diagonal part of the P-kernel of L, {q,Lq') = ¥{q)a{q\q'), cannot 
be negative. The smallest possible choice of S has as off-diagonal elements what is 
needed to compensate the negative values, and this leads to the minimal jump process, 
as described in Section 15.21 The diagonal part contains only what is needed to ensure 
that 5*1 = 0. For H of the form Hq + Hj, the role of S is again to compensate negative 
values off the diagonal, and the minimal process has velocities determined by Hq via 
(jn21) and jump rates determined by Hj via 

In any case, the backward generator of the minimal process is the one closest, in a 
suitable sense, to ()118p . This formula may thus be regarded as containing the essential 
structure of L, for the deterministic as well as for the jump part of the process. 

Another approach towards a general notion of minimal process may be to approxi- 
mate H by Hilbert-Schmidt operators Hn, with which are associated, according to the 
results of Sections 4.2.1 and 4.2.4 of [15j, minimal jump processes Qn, and take the limit 
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77, — >■ cxo of the processes Qn- This leads to a number of mathematical questions, such 
as under what conditions on H,"^, P, and if„ does a limiting process exist, and is it 
independent of the choice of the approximating sequence Hn- 

6 Remarks 
6.1 Symmetries 

Process additivity preserves symmetries, in the sense that the process generated by 
-^^^^ shares the symmetries respected by all of the building blocks =Sf This section 
elaborates on this statement, and the following ones: The minimal jump rates (j^^ and 
the minimal free generator ()32|) share the symmetries of the Hamiltonians with which 
they are associated. The "second quantization" algorithm preserves the symmetries 
respected by the one-particle process. 

Here are some desirable symmetries that may serve as examples: space translations, 
rotations and inversion, time translations and reversal, Galilean or Lorentz boosts, global 
change of phase \E' e*^\E', relabeling of particles, and gauge transformations. 

We focus first on symmetries that do not involve time in any way, such as rotation 
symmetry. In this symmetry group G acts on Q, so that to every g G G there 

corresponds a mapping (f^ : Q ^ Q. In addition, G acts on through a projective 
unitary (or anti-unitary) representation, so that to every g E G there corresponds a 
unitary (or anti- unitary) operator Ug. Then the theory is G-invariant if both the wave 
function dynamics and the process on Q are, i.e., if H is G-invariant, 

U;'HUg = H, (119) 

and 

^^(Qf)=gf^* (120) 

in distribution on path space. A necessary condition for p2()|l is that the "configuration 
observable" transforms like the configuration, in the sense that 

U^'P{-)Ug = ^lP{-), (121) 

where cp^ denotes the action of if on measures. Without p21|l . p2()j) would already fail 
at time t = 0, no matter what the generator is. Given ()12H) . p20|) is equivalent to the 
G-invariance of the generator: 

^^^Vr=if^«*. (122) 

Since (p^ is a linear operator, it follows immediately that the sum of G-invariant gen- 
erators is again G-invariant. The minimal jump process, when it exists, is G-invariant, 

^•^This may mean two things: changing the artificial labels given to identical particles, or exchanging 
two species of particles. 
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as follows from the fact that ip^a"^ {dq\(p^{q')) = a^'^'^ [dq\q') , which can be seen by in- 
specting the jump rate formula (j29j) . The minimal free generator (j35p satisfies 
by virtue of (jll9j) and (jl2H) . "Second quantization" provides G-actions on TQ^^^ and 
^ = TJif^'^^ from given actions on Q(^) and (HI^ . (fnT|) and (fT^ are inherited 

from their 1-particle versions. 

Time-translation invariance is particularly simple. Consider generators which 
do not depend on time except through their dependence on \E'. Then the same is true of 
^ _^(«) _ xhe same can be said of the "second quantized" generator, and, provided H is 
time-independent, of the minimal jump rates (j^^ and the minimal free generator ()35|) . 

Next we consider time reversal. It is represented on Jif by an anti-unitary operator 
T, i.e., an anti-linear operator such that (T$|T\E') is the conjugate of ($|\1'). We assume 
that the Hamiltonian is reversible, THT~^ = H. Then the reversibility of the theory 
means that 

Q":, = Qr (123) 

in distribution on path space, where the superscript should be understood as indicating 
the state vector at t = 0. The necessary condition analogous to ()121|) reads 

T-^P(-)T = P(-), (124) 

and given that, ()123|) is equivalent to the T-invariance of the generator: 

= ■^T'i ) or L-q, = Lx-^ , (125) 

where .if and L denote the forward and backward generator of the time-reversed process. 
L can be computed from L, for an equivariant Markov process, according to^'' 

y = L^f - iL^l)f (126) 

^"^To make this formula plausible, it may be helpful to note that the second term on the right hand 
side is just the correction needed to ensure that Ltl = 0, a necessary condition for being a backward 
generator. If P were stationary, the second term on the right hand side would vanish. 

Here is a derivation of H12()l) : Let {f,g) — J^^q f{q) g{q)P{dq) be the scalar product in i^(Q,P). It 
follows from the definition © of L that 

(5, Lf) = lim i E{g{Qo)fm - 5(Qo)/(Qo)) ■ 

Correspondingly, L is characterized (for / and g real) by 

(5,L/) = \imjE{g{Qo)fiQ^t)-9{Qo)f{Qo)) = 

= limjE{giQo)f{Q-t)~giQ-t)fiQ-t)) + 

+ lim i E(g(Q_0/(Q-0 - 9iQo)f{Qo)) = 

= (/, Lg) - I g{q) f{q) P(dg) 1* [Lg, /) - {L{gf), 1) = (5, L^) - (/g, LH) , 
96S 

which amounts to H126|) . 
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where ^ denotes the adjoint operator on L'^{Q, P), with P given by ()24|) . Since L is hnear 
in L, condition (jl25p is preserved when adding (forward or backward) generators; it is 
also preserved under "second quantization." For a pure jump process, (|125|) boils down 
to 

(^|P(dg')l*) = (^^""(dq'lq) (^|P(cig)|^) , (127) 

which is satisfied for the minimal jump rates, by inspection of ((221) • The minimal free 
generator (j32p changes sign when replacing \E' by T\E', which means the velocity changes 
sign, as it should under time reversal (see Section l^!T|) . 

Invariance under Galilean boosts is a more involved story, and as it is not considered 
as fundamental in physics anyway, we omit it here. Lorentz boosts are even trickier, since 
for more than just one particle, they even fail to map (simultaneous) configurations into 
(simultaneous) configurations. As a result, the problem of Lorentz invariance belongs 
in an altogether different league, which shall not be entered here. 

6.2 On the Notion of Reversibility 

It may appear, and it is in fact a widespread belief, that stochasticity is incompati- 
ble with time reversibility. We naturally view the past as fixed, and the future, in a 
stochastic theory, as free, determined only by innovations. Even Bell expressed such a 
belief P- 177]. However, from the proper perspective the confiict disappears, and this 
perspective is to consider the path space (of the universe) and the probability measure 
thereon. If t Qt is a. history of a universe governed by a Bell- type QFT, then its time 
reverse, 1 1— > Q-j, is again a possible path of this Bell-type QFT, though corresponding 
to a different initial state vector T\E' instead of \E', with T the time reversal operator as 
discussed in Section More than this, the distribution of the reversed path t i— > Q-t 
coincides with the probability measure on path space arising from T^.^^ 

It may also be helpful to think of how the situation appears when viewed from outside 
space-time: then the path Qt corresponds to the decoration of space-time with a pattern 
of world lines, and this pattern is random with respect to a probability measure on what 
corresponds to path space, namely the space of all possible decorations of space-time. 
Then time reversal is a mere reflection, and for a theory to be time reversible means 
the same as being invariant under this reflection: that we could have had as well the 
reflected probability measure, provided we had started with T\E' instead of \E'. 

To sum up, we would like to convey that the sense of reversibility for Markov pro- 
cesses indeed matches the sense of reversibility that one should expect from a physical 
theory. 

^^We can be more precise about the meaning of the measure on path space: as in Bohmian mechanics 
its role "is precisely to permit definition of the word 'typical'." jHl P- 129] Consequently, the 
meaning of the reversibility property of the measures we just mentioned is that the time reverse of a 
history that is typical with respect to ^, is typical with respect to T"^. 
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6.3 Heisenberg Picture 

In (1^ . we have applied the Schrodinger picture, according to which the state vector 
evolves while the operators remain fixed. Eq. and the reasoning following it can 
as well be translated to the Heisenberg picture where the state vector \E' is regarded as 
fixed and the operators Pt( • ) as evolving. Thus, we could equivalently write 

¥t{dq) = {^\Pt{dq)\^) 

instead of ()24|) . Similarly, Hq and Hj become time-dependent while their sum is con- 
stant. We often use an ambiguous notation like {'^\P{dq)\'^) and formula since 
the formulas are equally valid in both pictures (and, for that matter, in the interaction 
picture). 

Like the jump rate formula (j29p . the formula (j32j) for the minimal free generator is 
equally valid in the Heisenberg picture. 

We further remark that in the Heisenberg picture, the following nice equation holds 
for a pure jump process with minimal rates when P is a PVM: 

PT0h{Qt+at e dq,Qt G dq'} = {^\{Pt+,t{dq),Pt{dq')}\^)^ (128) 

for dq n dq' = 0, where { , } on the right hand side means the anti-commutator. The 
similarity to the one-time distribution formula 

Froh{Qt e dq} = {^\Ptidq)\^) 

is striking. Specifying the two-time distribution for infinitesimal time differences is a way 
of characterizing a Markov process, equivalent to specifying the (forward or backward) 
generator and the one-time distribution. Thus, for a PVM P (jl28|) provides another 
formula for the minimal jump rates ()29|). A similar formula for the process generated 
by the minimal free generator §^ is E{g{Qt)fiQt+dt)) = ^{"^{{gt, ft+dtW- 

6.4 Examples of Process Additivity 

Among different versions of Bohmian mechanics we find numerous examples of process 
additivity (and, remarkably, no example violating it): 

• The Hamiltonian for n noninteracting particles is the sum of the Hamiltonians 
for the individual particles, and it is easy to see that the vector field defining 
Bohmian mechanics for the ra-particle system is the sum of the vectors fields (each 
regarded as vectors fields on R^") for the particles. As already mentioned, sums 
of generators for deterministic processes amount to sums of the defining vector 
fields. 

Moreover, the vector field for each particle is essentially the Bohmian one-particle 
law. To point out that this is a nontrivial fact, we mention that this is not so for 
the alternative velocity formula (10.2) in |T2] considered by Deotto and Ghirardi, 
for which the velocity of the i-th particle differs from the one-particle law. So 
Bohmian mechanics of n particles can be viewed as built from n copies of the 
one-particle version, in fact by the "second quantization" algorithm of Section 
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• The vector field of Bolimian meclianics for a single spinless particle may also be 
seen as arising in this way. If a Hamiltonian H = —X"^ is the negative square 
of an (incompressible) vector field (regarded as a first-order differential operator) 
X = a{x)-V on (with V-a = ensuring formal self-adjointness of the square), 
then the simplest equivariant process associated with H is given by the velocity 
vector field 

2 ^ a ■ 

f = — Im a . 

h ^ 

The corresponding backward generator is L = |lm(^)X. Now = 
— J2a-^a^ °f 3 negative squares of vector fields Xa = -^d/dx^ corre- 

sponding to the individual degrees of freedom. The associated Bohm velocity is 
the sum of the velocities corresponding to the squares. So Bohmian mechanics in 
three dimensions can be viewed as built of 3 copies of the one-dimensional ver- 
sion. To point out that this is a nontrivial fact, we mention that this is not true, 
e.g., of the velocity formulas (10.1) and (10.2) in [12], which do not make sense in 
dimensions other than 3. 

• If we add an interaction potential V to — -yA, the Bohm velocity is the appropriate 
sum, since the operator V is associated with the trivial motion w = 0. 

• We may also include an external vector potential A{x, t) in the Schrodinger equa- 
tion, that is, replace -f A = -f by -f (V + 2|A(a;, t))^ = -f A-f (^|V• 
A + ^fA ■ V) + f The sum of the associated velocities, namely 

v|/*V\I^ 
film— — — + eA + 

equals the velocity one obtains directly, hlm^*{y + i|A)\E'/\E'* . 

• In the Bohm-Dirac theory (|31|). however, one can include an external gauge con- 
nection A^{x, t) in the Dirac equation without changing the velocity formula. That 
conforms with process additivity because the operator {'^^)~'^'^'^A^ = + a ■ A 
is associated (termwise) with = 0. 

• In the Dirac Hamiltonian H = —icha ■ V -|- Pmc^, the first term corresponds to 
the Bohm-Dirac velocity (j31|) . whereas the second term corresponds to f = 0; as 
a consequence, the Bohm-Dirac velocity does not depend on the mass. Moreover, 
the three components of the Bohm-Dirac velocity are each equivariant with respect 
to the corresponding derivative term in H. 

In addition, we point out cases of process additivity in the "second quantization" 
algorithm and minimal jump processes. 

The "second quantized" generator TJ^^^^ as constructed in Section provides an 
example of process additivity (or may be viewed as an application of process additivity): 



n=0 
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where the generators in the sum correspond to motions in the respective different sectors 
of Q. 

Suppose we regard the particles as ordered, Q = {Q^, . . . , Qj^). Then another case 
of process additivity becomes visible: 

where hi is the one-particle Hamiltonian acting on the i-th particle. Correspondingly, 

N 
i=l 

where =Sfj is equivariant with respect to hi. This applies not only to Bohmian me- 
chanics (as described earher in this section), but generally to the "second quantiza- 
tion" procedure as described in Section 14 .21 We also note that the "second quantiza- 
tion" algorithm presented in Section 14.21 preserves process additivity in the sense that 
r(^/') + ^2^')) = T{^i^^) + r(^2^')) while r(i7f ) + H^^^) = T{Hi^^) + T{Hi^^). 

We now turn to process additivity among minimal jump processes. 

A jump process generated by a sum need not be a minimal jump process even 
when its constituents are. But under certain conditions it is. Two such cases are the 
"direct sum" and "tensor product" processes constructed in Sections 14.2.11 and 14.2.21 

= 0„ ^(") with Q = Un 2^"^ and = 0„ and = ^t^l ® ■ • ■ ® ^1^1 

with Q = X ■ ■ ■ X and = ^ . 1 ® ■ • • ® H^^ ® ■ ■ ■ ® 1, with ^ = J2-^i where 
=Sfj acts nontrivially, in an obvious sense, only on Q*^*-* or on Q^*'. These are special cases 
of the general fact that minimality is compatible with additivity whenever the addends 
of the Hamiltonian correspond to different sorts of jumps. That can be most easily 
understood in the case of a PVM corresponding to an orthonormal basis {\q) '■ q & Q} 
of J^: suppose H = H1 + H2 and for every pair q, q' either {q\Hi\q') = or {q\H2\q') = 0. 
Then a = o"i + cr2. The corresponding condition in the POVM context is that the kernels 
of Hi and H2 have disjoint supports. When H is naturally given as a sum this condition 
would be expected to be satisfied. 

Finally, we remark that the minimal free generator ^ = as defined in fl35|) is 
additive in H. 

6.5 Second Quantization of a Minimal Jump Process 

We note that the "second quantization" of a minimal jump process associated with a 
PVM P*^^-*, as described in Section f4.2t is the minimal jump process associated with 
the second-quantized Hamiltonian; this is a consequence of the observation that S£i 
generates the minimal jump process for Hi in this case. This fact is probably physically 
irrelevant but it is mathematically nice. 
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6.6 Global Existence Question 

The rates at and velocities Vt, together with Pt, define the process Qt associated with 
H, P, and which can be constructed along the lines of Section 12.31 However, the 
rigorous existence of this process, like the global existence of solutions for an ordinary 
differential equation, is no trivial matter. See Section 4.3 of for a discussion of what 
must be controlled in order to establish the global existence of the process, and [20] for 
an example of such a global existence proof. 

6.7 POVM Versus PVM 

As we have already remarked in footnote El every POVM P is related to a PVM Pcxt, 
the Naimark extension, on a larger Hilbert space J^ext according to P( ■ ) = P+Pcxt( ■ )I 
with P+ the projection Jifext ^ and / the inclusion Jif J^xt- This fact allows 
a second perspective on P, and sometimes creates a certain ambiguity as to which 
process is the suitable one for a Bell-type QFT, as follows. At several places in this 
paper, we have described considerations leading to and methods for defining Markov 
processes, in particular minimal jump rates (j29|) and the minimal free generator (j32|) : 
these considerations and methods could be applied using either and Pext or and 
P. One would insist that the state vector must lie in Jif, the space of physical states, 
but even then one might arrive at different processes starting from P or Pext- To obtain 
a process from Pext requires, of course, that we have a Hamiltonian on ^xt, while H 
is defined on Jif; such a Hamiltonian, however, can easily be constructed from H by 
setting ifext = IHP+. 

In some cases, the Naimark extension does not lead to an ambiguity. This is the case 
for the jump rate formula since for e Jf, (^'|Pext(c?g)|^) = (^|P(rfg)|^) and 

(^|Pext('ig)-ffextPext('^g')l^) = | P((ig)-ff P(o?g') I • This fact suggests that, generally, 
the minimal process arising from i/ext and Pext is the same as the one arising from H 
and P. 

The situation is different, however, when H is defined on Jifext to begin with, and 
different from -ffext- This is the case with the free Dirac operator ho, defined as a 
differential operator on L^(M^,C^), which differs from P+/;,oP+. When we obtained in 
Section 12.91 the Bohm-Dirac motion (j31|) from the formula (j32p for the minimal free 
generator, we used and Pext- In contrast, the restriction of to the positive energy 
subspace, or equivalently P^hoP^, possesses a kernel; more precisely, it is a convolution 
operator S+ -k {hoS+)* in the notation of Section 13.31 and thus corresponds to jumps. 
The associated minimal process on presumably makes infinitely many jumps in every 
finite time interval, similar to the example of ^^1; Section 3.5. 

Thus, there are two processes to choose between, the Bohm-Dirac motion and the 
minimal process for P+hoPj^. Both are equivariant, and thus it is arguably impossible 
to decide empirically which one is right. In our example theory in Section ESI we chose 
the simpler, deterministic one. But we leave to future work the discussion of which is 
more likely relevant to physics, and why. 
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6.8 The Role of Field Operators 



The Bell-type QFTs with which we have been concerned in this paper are models de- 
scribing the behaviour of particles moving in physical 3-space, not of fields on 3-space. 
We have been concerned here mainly with a particle ontology, not a field ontology. This 
focus may be surprising at first: almost by definition, it would seem that QFT deals 
with fields, and not with particles. Consider only the occurrence (and prominence) of 
field operators in QFT! 

But there is less to this than might be expected. The field operators do not function 
as observables in QFT. It is far from clear how to actually "observe" them, and even 
if this could somehow, in some sense, be done, it is important to bear in mind that 
the standard predictions of QFT are grounded in the particle representation, not the 
field representation: Experiments in high energy physics are scattering experiments, in 
which what is observed is the asymptotic motion of the outgoing particles. Moreover, for 
Fermi fields — the matter fields — the field as a whole (at a given time) could not possibly 
be observable, since Fermi fields anti-commute, rather than commute, at space-like 
separation. One should be careful here not to be taken in by the attitude widespread 
in quantum theory of intuitively regarding the operators as "quantities," as if they 
represented something out there in reality; see P for a critique of this attitude. 

So let us focus on the role of the field operators in QFT. This seems to be to relate 
abstract Hilbert space to space-time: the field operators are attached to space-time 
points, unlike the quantum states \E', which are usually regarded not as functions but as 
abstract vectors. In orthodox quantum field theory the field operators are an effective 
device for the specification of Hamiltonians having good space-time properties. For our 
purposes here, what is critical is the connection between field operators and POVMs. 

Throughout this paper, the connection between Hilbert space and the particle posi- 
tions in physical space has been made through the POVM P, and through it alone. We 
now wish to emphasize that the field operators are closely related to P, and indeed that 
field operators are just what is needed for efficiently defining a POVM P on r(R'^). 

This connection is made through number operators N{R), R C M^. These define 
a number- operator-valued measure (NOVM) A'"( • ) on R^, an "unnormalized POVM" 
(iV(M^) is usually not the identity operator and N{R) is usually an unbounded 
positive operator) for which the values N{R) commute and are number operators: 
spectrum(A^(i?)) C {0, 1,2,3, . . .}. (The basic difference, then, between a NOVM and 
a PVM is that the spectrum of the positive operators is {0, 1, 2,3,.. .} rather than just 



There is an obvious one-to-one relation between NOVMs N{ ■ ) on and PVMs P 
on r(]R^), given by 



where nji{q) = #(g fl R) is the number function on r(R^) for the region R. Since (I129|) 
is the spectral decomposition of the commuting family N{R), this correspondence is 
one-to-one. (Note that the joint spectrum of the commuting family N{R) is the set of 



{0,1}.) 




(129) 



ger(M3) 
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nonnegative-integer-valued measures on M^, one of the definitions of r(M^) given in 
Section EHl) 

The moral is that a NOVM on is just a different way of speaking about a PVM 
P on Q = r(M^). All other POVMs arise from PVMs by restriction to a subspace 
(Naimark's theorem ^I]). An easy way to obtain a NOVM starts with setting 



for suitable operators (l){x). An easy way to ensure that the N{R) commute is to require 
that the operators (pix) commute or anti-commute with each other and the adjoints 
4>*{x') for x' 7^ X. An easy way to ensure that the N{R) have nonnegative integer 
eigenvalues is to require that 



where [ , ]± is the (anti-)commutator, and that there is a cyclic vacuum state |0) G 
for which (/)(£c)|0) = 0. The relations 1)1311) are of course just the usual canonical 
(anti-)commutation relations that field operators are required to satisfy. 

Moreover, in gauge theories the connection between matter field and the NOVM 
is perhaps even more compelling. Consider a gauge theory with internal state space 
V, equipped with the inner product ((-I-)). Then, given x G M'^, the matter field 
(f){x) should formally be regarded as a linear functional V —* 0{Jif), ^ (j)^{x), from 
the internal state space to operators on Jif, with (as) = {(j)^{x))* a linear function 
V* 0{J^) on the dual of V. (HSH) then becomes [0^(a;), 0;.(a;')] = S{x - x') {{7]\Q. 
Thus the simplest gauge-invariant object associated with is the NOVM (I13(jp . with 
the integrand understood as the contraction of the tensor V x V* ^ 0{M'), {i^rf) i— >• 



Hence, not only does the notion of particle not conflict with the prominence of field 
operators (see Sections 13.11 and 13.31 for exphcit examples), but field operators have a 
natural place in a theory whose ultimate goal it is to govern the motion of particles. 
One of their important roles is to define the POVM P that relates Hilbert space to 
configuration space. Quantum theory of fields or quantum theory of particles? A theory 
of particle motion exploiting field operators! 

7 Conclusions 

The essential point of this paper is that there is a direct and natural way of under- 
standing QFT as a theory about moving particles, an idea pioneered, in the realm of 
nonrelativistic quantum mechanics, by de Broglie and Bohm. We leave open, however, 
three considerable gaps: the question of the process associated with the Klein-Gordon 
operator, the problem of removing cut-offs, and the issue of Lorentz invariance. 
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